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Abstract. This is the second in a series of papers where we combine the 
classical approach to exponential Diophantine equations (linear forms in loga- 
rithms, Thue equations, etc.) with a modular approach based on some of the 
ideas of the proof of Fermat's Last Theorem. In this paper we give a general 
and powerful lower bound for linear forms in three logarithms. We use this 
lower bound, together with a combination of classical, elementary and sub- 
stantially improved modular methods to solve completely the Lebesgue-Nagell 
equation 

+ D = y" , X, y integers, n > 3, 
for D in the range I < D < 100. 



1. Introduction 

Arguably, the two most celebrated achievements of the 20th century in the field of 
Diophantine equations have been Baker's theory of linear forms in logarithms, and 
Wiles' proof of Fermat's Last Theorem. We call Baker's approach to Diophantine 
equations the 'classical approach'. The proof of Fermat's Last Theorem is based on 
what we term the 'modular approach'. The proponents of the classical approach 
are too many to mention; the modular approach is still in its infancy, but among 
the early contributers let us just mention Frey, Serre, Ribet, Darmon, Merel, Kraus, 
Bennett, Skinner, Ivorra, etc. 

The motivation for our series of papers, of which this is the second, is that neither 
approach (on its own, and as it stands at the moment) is powerful enough to resolve 
unconditionally many of the outstanding exponential Diophantine equations. Our 
thesis is that one should, where possible, attack exponential Diophantine equations 
by a combination of the classical and modular approaches. The precise aims of this 
series were formulated in our first paper [H] as follows: 

(I) To present theoretical improvements to various aspects of the classical ap- 
proach. 

(II) To show how local information obtained through the modular approach 
can be used to reduce the size of the bounds, both for exponents and for 
variables, of solutions to exponential Diophantine equations. 
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(III) To show how local information obtained through the modular approach can 
be pieced together to provide a proof that there are no missing solutions 
less than the bounds obtained in (I), (II). 

(IV) To solve various famous and hitherto outstanding exponential Diophantine 
equations. 

In 1^ we gave a lower bound for linear forms in three logarithms, and used a 
combination of classical and modular methods to determine all the perfect powers 
in the Fibonacci and Lucas sequences. In this paper, we give a new lower bound 
for linear forms in three logarithms that is more general and powerful than the one 
given in the previous paper. We are also concerned with the following exponential 
Diophantine equation, which we call the Lebesgue-Nagell equation: 

(1) + D = y", X, y integers, n> 3. 

Here, D denotes a non-zero integer. The reason for the name Lebesgue-Nagell is 
given in Section El together with some historical remarks. But for now we men- 
tion that the equation had previously been solved for 81 values of D in the range 
1 < D < 100, using elementary, classical and modular methods; the remaining 
values are clearly beyond these methods as they stand. In this paper we apply 
our lower bound for linear forms in three logarithms, together with a combination 
of elementary, classical, and substantially improved modular methods to prove the 
following Theorem. 

Theorem 1. All solutions to equation with D in the range 

(2) 1<D< 100 

are given in the Tables at the end. In particular, the only integer solutions {x, y, n) 
to the equation 

x'^ + 7 = y", n > 3, 
satisfy \x\ = 1, 3, 5, 11, 181. 

We choose to give a complete proof of Theorem 1, rather than treating the 19 
remaning values of D in the range ^ . 

It is noted that the solutions for even n can be deduced quickly, for then D is 
expressible as a difference of squares. It is therefore sufficient to solve the equation 

(3) x^ + D = yP, X, y integers, p > 3 is prime; 

the solutions to Q can then be recovered from the solutions to Q . 

We give three modular methods for attacking Two are refinements of known 
methods, and a third that is completely new. Using a computer program based 
on these modular methods, we can show - for any D in the above range - that 
the exponent p is large (showing that p > 10^ is quite practical). Our modular 
approach also yields the following rather surprising result: either each prime factor 
of y divides 2D, or y > (y/p — 1)^. We are then able to deduce not only that p is 
large, but also that y is large. This information helps to reduce the size of the upper 
bound on p obtained from the lower bound for the linear forms in three logarithms, 
making the computation much more practical. Our total computer time for the 
computations in this paper is roughly 206 days on various workstations (the precise 
details are given in due course). 
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Using our approach should make it possible to solve for any reasonable D 
that is not of the form D ~ —a} ± 1; if is of this form then the equation has 
a solution (a;,y) = (a, ±1) for all odd values of the exponent n, and the modular 
methods we explain later are not very successful in this situation. To deal with 
this case requires further considerations which we leave for another paper. Notice, 
however, that we solve the case D = 1. 

We would like to warmly thank Mihai Cipu for pointing our many imperfections 
in a previous version of this paper, and Guillaume Hanrot for help with solving 
Thue equations. 

2. On the History of the Lebesgue-Nagell Equation 

Equation JQ) has a long and glorious history, and there are literally hundreds 
(if not thousands) of papers devoted to special cases of this equation. Most of 
these are concerned with equation |^ either for special values of n or special values 
of y. For example, for 13 = 2, n = 3, Fermat asserted that he had shown that 
the only solutions are given by x = 5, ?; = 3; a proof was given by Euler 
Equation ^ with n = 3 is the intensively studied Mordell equation (see fo^' ^ 
modern approach). 

Another notable special case is the generalized Ramanujan-Nagell equation 

(4) a;2 + i:> = fc", 

where D and k are given integers. This is an extension of the Ramanujan-Nagell 
equation + 7 = 2", proposed by Ramanujan [331 in 1913 and first solved by 
Nagell in 1948 (see also the collected papers of Nagell |33J). This equation has 
exactly five solutions with a; > 1 and is in this respect singular: indeed, Bugeaud and 
Shorey ^01 established that equation with D positive and fc a prime number 
not dividing D has at most two solutions in positive integers n, except for 
{D,k) = (7,2). They also list all the pairs {D,k) as above for which equation Q 
has exactly two solutions. We direct the reader to for further results and 
references. 

Returning to equation the first result for general y, n seems to be the proof in 
1850 by V. A. Lebesgue ^22, that there are no non-trivial solutions for D = 1. The 
next cases to be solved were D = 3, 5 by Nagell [33] in 1923. It is for this reason 
that we call equation the Lebesgue-Nagell equation. The case with D = — 1 is 
particularly noteworthy: a solution was sought for many years as a special case of 
the Catalan conjecture. This case was finally settled by Chao Ko [21] in 1965. 

The history of the Lebesgue-Nagell equation is meticulously documented in an 
important article by Cohn jT3] . and so we are saved the trouble of compiling an 
exhaustive survey. In particular, Cohn refines the earlier elementary approaches 
of various authors and completes the solution for 77 values of D in the range 
1 < D < 100. The solution for the cases D = 74, 86 is given by Mignotte and de 
Weger [221 ■ Bennett and Skinner [HJ Proposition 8.5] apply the modular approach 
to solve D = 55, 95. The 19 remaining values 

(5) 7, 15, 18, 23, 25, 28, 31, 39, 45, 47, 60, 63, 71, 72, 79, 87, 92, 99, 100, 

are clearly beyond the scope of Cohn's elementary method, though Cohn's method 
can still give non-trivial information even in these cases, and is revisited in Section]^ 
Moreover, as far as we can see, the modular method used by Bennett and Skinner 
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(which is what we later on caU Method I) is not capable of handling these values 
on it own, even though it still gives useful information in most cases. 

Cohn, in the same paper, also makes a challenge of proving that the only solutions 
to the equation 

x^ + 7 = y" 

have l^l = 1, 3, 5, 11, 181. This challenge is taken up by Siksek and Cremona 
021 who use the modular approach to show that there are no further solutions 
for n < 10^, nor for composite n. They also suggest that an improvement to lower 
bounds in linear forms in three logarithms may finally settle the problem. With the 
benefit of hindsight, we know that they were almost - though not entirely - correct. 
The substantial improvement to lower bounds in linear forms in three logarithms 
given here, was certainly needed. However, for this lower bound to be even more 
effective, a further insight obtained from the modular approach was also needed: 
namely that y is large as indicated in the introduction. 

3. Reduction to Thue Equations 

Our main methods for attacking equation ^ are linear forms in logarithms (to 
bound the exponent p) and the modular approach, though for some small values 
of p it is necessary to reduce the equation to a family of Thue equations. The 
method for reducing equation Q to Thue equations is well-known. We do however 
feel compelled to give a succinct recipe for this, in order to set up notation that is 
needed later. 

It is appropriate to point out that there are other approaches that could be used 
to solve equation (j^J for small p. For p = 3 we can view the problem as that of 
finding integral points on elliptic curve, a problem that is aptly dealt with in the 
literature (see 01] and ^7]). For p > 5, the equation + D = yP defines a curve of 
genus > 2; one can sometimes determine all rational points on this curve using the 
method of Chabauty JTj, though this would require computing the Mordell-Weil 
group of the Jacobian as well (see [HH], [HHI, ESI, EH and |1H|). 

We do not assume in this section that D is necessarily in the range (O, merely 
that — Z) is not a square. We write (here and throughout the paper) 

D — D\D2, Di, D2 are integers, D2 square-free. 

Let £ = Q{y/—D2), and O be its ring of integers. Throughout the present paper, 
we denote the conjugate of an element a (resp. of an ideal a) by a (resp. by a). 

Let pi, . . . , pr be the prime ideals of O dividing 2D. Let A be the set of integral 
ideals a of O such that 

• a = pi^ • • • p^-vwith < a, < p, 

• (a, a) I 2Di^^D^, 

• the ideal ad is a perfect p-th power. 

If {x, y) is a solution to equation 0, then one effortlessly sees that 

{x + Di^^D2)0 = abP 

for some a £ A and some integral ideal b. 

Now let bi, . . . , b/i be integral ideals forming a complete set of representatives 
for the ideal class group of O. Thus bbi is a principal ideal for some i, and so 
bbi = (3'0 for some (3' € O. The fractional ideal ab~^ is easily seen to be also 
principal. The ideal b is unknown, but the ideals, a,bi, . . . ,bh are known. We may 
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certainly determine which of the fractional ideals ab^^ are principal. Let T' be a 
set containing one generator 7' for every principal ideal of the form ab~^ {a €z A 
and 1 < i < h). It is noted that the elements of T' are not necessarily integral, but 
we know that if {x, y) is a solution to equation Q then 



for some 7' G F', and some (3' e O. Finally, define F as follows: 

f F' if D2 >0,D2^ 3, or if L>2 = 3 and p ^ 3, 

F = i F' U Cr' U C^^r' if Da = p = 3, where C = (1 + V^)/2, 

[Uje-^F' if D2 < 0, where j ranges over -{p - l)/2, . . . , (p - l)/2, 

where if D2 < (and so C is real) we write e for the fundamental unit. 
We quickly deduce the following. 

Proposition 3.1. With notation as above, if {x,y) is a solution to equation ^ 
then there exists 7 € F and f3 d O such that 



x + Di^^^-fl3P. 
Thus if we let 1, oj be an integral basis for O then for some 7 € F, 

X = i (7([/ + VujY + 7(L/ + VluY) 

for some integral solution (U, V) to the Thue equation 

^ {j{U + Vlj)p -j{U + Vujy) = Di . 



3.1. Results I. If q is a prime we denote by : Z Z>o U {00} the normalized 
g-adic valuation. 

We now eliminate all cases where it is inconvenient to carry out level-lowering. 

Lemma 3.2. Suppose D is in our range 0j. Suppose (x, y,p) is a solution to equa- 
tion 0^ that is missing from our Tables at the end. Then p satisfies the following 
conditions: 

r p>7, 

(6) \ P > ^q{D) + 1 /c" 0,11 primes q, 

\ p> V2{D) + 7 if V2{D) is even. 

Proof. It is clear that for any particular D there are only a handful of primes p 
violating any of these conditions. We wrote a pari/gp 1 program that solved all 
the equations © for p violating © : the program first reduces each such equation to 
a family of Thue equations as in Proposition l3 . II above. These are then solved using 
the in-built pari/gp function for solving Thue equations (this is an implementation 
of the method of Bilu and Hanrot [S]). 

It is perhaps worthwhile to record here two tricks that helped us in this step. 
First, in writing down the set F appearing in Proposition 13.11 we needed a set of 
integral ideals bi, . . . , representing the ideal class group of the quadratic field £. 
Both pari/gp and MAGMA have in-built functions that amount to homomorphisms 
from the ideal class group as an abstract group, to the set of fractional ideals, and 
these can be used to construct the required set bi, . . . , b;,. We have found however 
that we get much simpler Thue equations if we search for the smallest prime ideal 
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representing each non-trivial ideal class, and of course taking lO to represent the 
trivial ideal class. 

To introduce the second trick, we recall that when one is faced with a Thue 
equation 

aoUP + aiUP-W + ■■■ + apV'P = b 

it is usual to multiply throughout by Aq ^ and make the substitution U' — a^U , 
thus obtaining a monic polynomial on the left-hand side. When oq is large, this 
greatly complicates the equation. The second trick is to first search for a unimodular 
substitution which makes the leading coefficient oq small. 

After optimizing our program, we were able to complete the proof in about 22 
minutes on a 1050 MHz UhraSPARC III computer. □ 

4. Removing Common Factors 

It is desirable when applying the modular approach to equation to remove 
the possible common factors of the three terms in the equation. This desire leads 
to a subdivision of cases according to the possible common factors, as seen in the 
following elementary Lemma. Here and elsewhere, for a non-zero integer o, the 
product of the distinct prime divisors of a is called the radical of a, and denoted 
by rad(a). 

Lemma 4.1. Suppose that (x, y,p) is a solution to equation 0) such that y ^ and 
p satisfies the condition Then there are integers di, d2 such that the following 
conditions are satisfied: 

(i) di > 0, 

(ii) D = dld2, 

(iii) gcd(di,d2) = 1, 

(iv) for all odd primes q \ di we have ^ — 1, 

(v) if2\di then d2 = 7 (mod 8). 
Moreover there are integers s, t such that 

X = dit, y = rad((ii)s, 

where rad(di) denotes the radical of di, and 

(7) t^ + d2 = esP, gcd(t, d2) = 1, s 7^ 0, 
where 

(8) e= ]J qP-2".(rfi), 

q prime 
q\di 

and rad(e) = rad((ii). 

Proof. Suppose {x,y,p) is a solution to equation Q such that y ^ and condi- 
tion lO is satisfied. It is straightforward to see that condition ^ forces gcd{x^, D) 
to be a square, say d\ with di > 0. We can therefore write x = dit and D = dld2 
for some integers t,d2. Moreover, since 

dj = gcd(a;^ D) = gcdidjt^, dld2) = dl gcd{t\d2), 

we see that gcd(t, ^2) — 1. Removing the common factors from x^ + D — yP we 
obtain t^ + d2 = esP where e is given by 0. The integrality of e follows from the 
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condition ©, and so does the equality of the radicals rad(e) — rad((ii). Note that 
(iii) follows from this equality of the radicals and the fact that t, c?2 are coprime. 
We have thus proven (i) , (ii) , (iii) and it is now easy to deduce (iv) and (v) . Finally, 
the condition s 7^ follows from the condition y ^ 0. □ 

Definition. Suppose D is a non-zero integer and {x,y,p) is a solution to equa- 
tion 0) with y ^ and p satisfying Let c?i, ^2 be as in the above Lemma and 
its proof (thus di > and gcd(x, D) — d\ and c?2 — D /d\). We call the pair {di, ^2) 
the signature of the solution {x,y,p). We call the pair {t,s) the simplification of 
{x,y) (or {t,s,p) the simplification of {x,y,p)). 

In this terminology, Lemma |4 . 1 1 associates to any D a finite set of possible signa- 
tures (c?i, c?2) for the solutions [x, y,p) of equation Q satisfying © and y ^ 0- To 
solve (PJ it is sufficient to solve it under the assumption that the solution's signature 
is (di,c?2) for each possible signature. 

Example 1. For example, if 13 = 25, there are two possible signatures satisfying 
the conditions of Lemma l4.ll these are (di,d2) — (1,25) or (5,1). If (c?i,o?2) = 
(1, 25), then x ^ t, y ^ s and we must solve the equation 

t^ + 25 = sP, 5\t. 

However, if (c?i, c?2) = (5, 1), then x — 5t, y — 5s, and we must solve the equation 

t^ + 1 = bP-'^sP. 

In either case it is noted that the three terms of the resulting equation are relatively 
coprime, which is important when we come to apply the modular approach. 



5. A Simplification of Cohn 



We will soon apply our modular machinery to equations Q with D in the 
range Q. Before doing this it is helpful to introduce a simplification due to Cohn 
that will drastically reduce the amount of computation needed later. All the argu- 
ments presented in this Section are found in Cohn's papers j ■ Cohn however 
assumes that D ^ 7 (mod 8), and the result that we state below is not formulated 
explicitly that way in these papers. 

Proposition 5.1. Let D — D\D2 where D2 is square-free and D2 > 0. Suppose 
that {x,y,p) is a solution to equation 0) with p satisfying J^, and let (c?i,c?2) be 
the signature of this solution. Then either 

(i) di > 1, 

(ii) or D EE 7 (mod 8) and 2 | y, 

(iii) or p divides the class number h of the quadratic field 



-D2), 



(iv) or y — a -\- D2b for some integers a, b such that 
b\D,, b^±Di, p\{Dl-b^), 
and a is a solution of the equation 

[u + b^^2Y -{u- b^^2r 



1 



(v) or D 



2V^ 
1. {x,y) 



(0,1), 
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(vi) or D2 = S (mod 4) and y — {a? + £'2&^)/4 for some odd integers a, b such 
that 

b\Di, p\{4Dl~b^), 
and a is a solution of the equation 
1 



{U + b^/^)P -{U- b^^2 



2PDi. 



Proof. We only give a brief sketch. Suppose that (i), (ii), (iii) are false. Then 
(x + Di-y/— -D2) = for some a in the ring of integers of Q(-\/— D2). There are two 
possibilities. The first is that a = a + by/—D2 for some integers a, b. By equating 
the imaginary parts we deduce all of (iv) ii b ^ ±Di. Thus suppose that b = ±Di. 
Letting (3 = a — b\/~D2 we see that 

a- 13 

If a/ (3 is not a root of unity, then the left-hand side is the p-th term of a Lucas 
sequence (with p > 7) and a deep Theorem of Bilu, Hanrot and Voutier [Sj on 
primitive divisors of Lucas and Lehmer sequences immediately gives a contradiction. 
Thus a/ (3 is a root of unity and so equal to ±1, ±i, or (±1 ± -\/— 3)/2. Each case 
turns out to be impossible, except for a = —j3 which together with b = ±Di implies 
(v). 

The second possibility for a is that a — {a + b\/—D2)/2 with a, b odd integers 
(and —D2 = 1 (mod 4)). Now (vi) follows quickly by equating the imaginary parts 

5.1. Results II. 

Corollary 5.2. Suppose D belongs to our range ^ and (x,y,p) is a solution to 
equation with p satisfying the condition If the solution (x,y,p) is missing 
from our Tables, then either D = 7 (mod 8) and 2 \ y, or di > 1 where {di,d2) is 
the signature of the solution. 

Proof. We apply Proposition 15. II Using a short MAGMA program we listed all so- 
lutions arising from possibilities (iv)-(vi) of that Proposition with \ < D < 100. 
The only ones found in our range are {x,y,p) = (0, l,p) for D — I and (x,y,p) = 
(±8, 2, 7) for D = 64 and these are certainly in the Tables. 

To prove the Corollary we merely have to take care of possibility (iii) of the 
Proposition. For 1 < D < 100, and primes p satisfying ©, the only case when p 
could possibly divide the class number of Q(\/— D2) is p = 7 and D = 71 (in which 
case h = 7). We solved the equation x'^ + 71 = y"^ by reducing to Thue equations as 
in Sectional It took pari/gp about 30 minutes to solve these Thue equations, and 
we obtained that the only solutions are {x,y) = (±46,3), again in our Tables. □ 

6. Level-Lowering 

In this section we apply the modular approach to equation Q under suitable, 
but mild, hypotheses. Ordinarily, one would have to construct a Frey curve or 
curves associated to our equation, show that the Galois representation is irreducible 
(under suitable hypotheses) using the results of Mazur and others jSO], and modular 
by the work of Wiles and others jSJi SHj; |H]i and finally apply Ribet's level- 
lowering Theorem Fortunately we are saved much trouble by the excellent 
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Table 1. Frey Curves with di, 6.2 odd. 



Case 


Condition on d2 


Condition on t 


Frey Curve Et 


L 


(a) 


d2 = l (mod 4) 




=X'^ + 2tX'^ - d2X 




(b) 


d2 = 3 (mod 8) 




y2 ^ x'A + 2tX'^ + + d2)X 




(c) 


d2 = 7 (mod 8) 


t even 


Y'^ ^ X'^ + 2tX'^ + (i^ + d2)X 




(d) 


d2 = 7 (mod 8) 


t = l (mod 4) 


Y"^ + XY - X^ + X"^ + ^ 


2 



Table 2. Frey Curves with d\ even, d2 odd. 



Case 


Conditions on s,p 


Frey Curve Et 


L 


(e) 


i = 1 (mod 4) 


y2 _^ xy = x3 + (*^) ^2 + 


^ 64 ; 




1 



paper of Bennett and Skinner which does aU of this for equations of the form 
Ax^ + By" = Cz^; it is noted that equation is indeed of this form. 

Let D be a non-zero integer. We shall apply the modular approach to the 
Diophantine equation 

(9) + £> = y^, \D, 0, and p > 3 is prime 
or the equivalent equation for the simplification (s,t) 

(10) t^+d2^esP, t^±l, gcd(t,rf2) = 1, s^O 

under the additional assumption that p satisfies ^ . The assumptions made about 
s, t in (|10|l are there to ensure the non-singularity of the Frey curves, and the 
absence of complex multiplication when we come to apply the modular approach 
later on. Before going on we note the following Lemma which in effect says that 
there is no harm in making these additional assumptions for D in our range 

Lemma 6.1. There are no solutions to the equation ^ for D in the range 0) 
with y — 0, or \ D, except those listed in the Tables at the end. 

Proof. Clearly y ^ 0. We produced our list of solutions with | D using a short 
MAGMA program. □ 

Lemma [4.11 associates to each equation of the form (jSJ finitely many signatures 
((ii,c?2) satisfying conditions (i)~-(v), and corresponding equations (0). Following 
Bennett and Skinner 3' we associate a Frey curve Et to any potential solution of 
equation H1Q(I according to Tables ^ |21 El 

The three tables are divided into cases (a)-(l). We know that di, d2 are coprime, 
and hence at most one of them is even. The possibility that di , d2 are both odd is 
dealt with in Tabled In cases (a), (b), a simple modulo 8 argument convinces us 
that t is odd. However for cases (c) and (d) - where di is odd and ^2 = 7 (mod 8) - 
the integer t can be either odd or even and we assign different Frey curves for each 
possibility. When t is odd (case (d)) we add the assumption that t = 1 (mod 4). 
This can be achieved by interchanging t with —t if necessary. 

Table deals with the possibility of even di, and Table 13 with the possibility of 
even d2. In both these cases t is necessarily odd, and the congruence condition on 
t can again be achieved by interchanging t with —t if necessary. 
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Table 3. Frey Curves with di odd, ^2 even. 



Case 


Condition on d2 


Condition on t 


Frey Curve Et 


L 


(f) 


V2{d2) = 1 




Y'^ = X'-* + 2tX'^ - d2X 


2« 


(g) 


^2 = 4 (mod 16) 


t = l (mod 4) 


Y2 = + tx^ - fx 


2 


(h) 


d2 = 12 (mod 16) 


t = 3 (mod 4) 


Y^ =X^ + tX^- fx 


22 


(i) 


"2(^2) = 3 


t = l (mod 4) 


r2 ^ x^ + tx^ ~ fx 


24 


(j) 


V2{d2) =4,5 


t = I (mod 4) 


r2 ^ ^3 + tx"^ - fx 


22 


(k) 


V2{d2) = 6 


t = l (mod 4) 


y2 + = x3 + (*^) x2 - ^x 


2-1 


(1) 


V2id2) > 7 


t = l (mod 4) 


Y"^ + XY = X^ + (i^) X2 - ^x 


1 



Proposition 6.2. Suppose D, di, d2 are non-zero integers that satisfy (i)-(v) of 
Lemma Suppose also p is a prime number satisfying the condition and 

let e he as defined in Suppose that (t, s) is a solution of equation and 

satisfying the supplementary condition (if any) on t in Tahles^\^\^ Let Et and L 
he as in these tables, and write Pp{Et) be the Galois representation on the p -torsion 
of Et- Then the representation pp(Et) arises from a cuspidal newform of weight 2 
and level N — Lva,d{D). 

Proof. The paper of Bennett and Skinner 3 gives an exhaustive recipe for Frey 
curves and level- lowering for equations of the form Ax"" + By"" = Cz^ under the as- 
sumption that the three terms in the equation are coprime. After a little relabeling, 
their results apply to our equation (|10|l and the Lemma follows from Sections 2, 3 
of their paper. It is here that we need the assumptions t 7^ ±1 and s 7^ made 
in (Cni) □ 

It is convenient to indulge in the following abuse of language. 

Definition, //(i, s,p) is a solution to equation J^l(]\) and if the representation pp{Et) 
arises from a cuspidal newform f , then we say that solution (t, s,p) arises from the 
newform f (via the Frey curve Et), or that the newform f gives rise to the solution 
{t,s,p). Lf (t,s,p) is the simplification of{x,y,p) then we say that {x,y,p) arises 
from the newform f. 

If the newform f is rational, and so corresponds to an elliptic curve E, then we 
also say that the solution {t,s,p) (or {x,y,p)) arises from E. 

6.1. A Summary. It may be helpful for the reader to summarize what we have 
done and where we are going. Given a non-zero integer D we would like to solve 
equation @. We can certainly write down all solutions with y = or with x'^ \ D. 
We can also solve (at least in principle) all cases where p violates condition JH)) by 
reducing to Thue equations as in Section |21 We can thus reduce to equation @ 
and assume that p satisfies condition ((HJ. 

Next, we can write down a list of signatures (^1,^2) satisfying conditions (i)-(v) 
of Lemma Em We reduce the solution of equation © to solving for each signature 
{di,d2) the equation (^UJ. Now we associate to the signature {di,d2) one or more 
Frey curves Et and levels L, so that any solution to (|10|l arises from some newform 
/ at level L via the Frey curve Et- 

Finally (and this is to come) we must show how to solve 1)10(1 under the assump- 
tion that the solution arises from a newform / via a Frey curve Et . If we can do 
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this for each ncwform / and Frey curve Et then we will have completed the solution 



As we shall see, the assumption that a solution arises from a particular newform 
is a very strong one, for it imposes congruence conditions on t modulo all but 
finitely many primes /. 

6.2. Congruences. For an elliptic curve E we write jJi?(F/) for the number of 
points on E over the finite field F;, and let ai{E) = Z + 1 — tti?(F;). 

Lemma 6.3. With notation as above, suppose that the Galois representation pp{Et) 
arises from a cuspidal newform with Fourier expansion around infinity 



of level N (given by Provosition 1 6. S]) and defined over a number field K/Q. Then 
there is a place ^ of K above p such that for every prime I \ 2pD we have 

ai{Et) = ci (mod *P) ift'^ + d2^0 (mod I) (or equivalently I \ s), 
I + 1 = ±ci (mod *P) ift'^ + d2 = (mod I) (or equivalently I \ s). 

Proof. The Lemma is standard (see jJDl, page 196], O page 7], [221 Proposition 
5.4], etc.). The conditions I \ 2D and I \ s together imply that I is a prime of good 
reduction for Et , whereas the conditions I \ 2D and I \ s imply that I is a prime of 
multiplicative reduction. □ 

When the newform / is rational, there is an elliptic curve E defined over Q 
whose conductor is equal to the level of the newform / such that a; (E) — ci for all 
primes I. In this case we can be a little more precise than in Lemma 16.31 thanks to 
a result of Kraus and Oesterle. 

Lemma 6.4. With notation as above, suppose that the Galois representation pp[Et) 
arises from a rational cuspidal newform f corresponding to an elliptic curve E/Q. 
Then for all primes I \ 2D we have 

ai{Et) = ai{E) (mod p) if t'^ + d2 ^ (mod I) (or equivalently I \ s), 
I + 1 = ztai{E) (mod p) if t^ + d2 = (mod I) (or equivalently I \ s). 

Proof. This Lemma does appear to be a special case of Lemma 16.31 however we do 
allow in this Lemma the case I — p which was excluded before. In fact Lemma 16.31 
together with a result of Kraus and Oesterle f35' Proposition 3] implies that the 
representations Pp{Et) and Pp{E) are semi-simply isomorphic. In this case the 
result of Kraus and Oesterle also tells us that ai{Et) = ai{E) (mod p) if the prime 
I is a prime of good reduction for both curves, and ai{Et)ai{E) = Z + 1 (mod p) if I 
is a prime of good reduction for one of them and a prime of multiplicative reduction 
for the other. Now since I \ 2D we see that / does not divide the conductor N oi E 
(which is also the level of the newform / as given by Proposition 16. 2|) . If Z ] s then 
I is a prime of multiplicative reduction for Et and then ai(Et) = ±1. The Lemma 
follows. □ 



of our equation lpj)l. 



(11) 




n>2 



7. Eliminating Exponents: Method I 



We now focus on equations of the form 1)1 0|l where, as always, p satisfies 
Proposition 16.21 tells us that if {t,s,p) is a solution to (|l()|l . then it arises from a 
newform of a certain level (or levels) and all these can be determined. Let us say 
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that these newforms are /i, . . . , /„. Then to solve equation (|10|l it is sufficient to 
solve it, for each i, under the assumption that the solution arises from the newform 
fi. We give three methods for attacking equation (|10|) under the assumption that 
the solution arises from a particular newform /. 

If successful, the first method will prove that the equation l|10|) has no solutions 
except possibly for finitely many exponents p and these are determined by the 
method. This method is actually quite standard. We believe that the basic idea is 
originally due to Serre jJU] pages 203-204] . It is also found in Bennett and Skinner 
[3 Proposition 4.3]. We shall however give a more careful version than is found in 
the literature, thereby maximizing the probability of success. 

Proposition 7.1. (Method I) Let D, di, d2 be a triple of integers satisfying (i)- 
(v) of Lemma \4.1\ Let f be a newform with Fourier expansion as in ^TJj having 
coefficients in the ring of integers of a number field K , and let Mk/q denote the 

q) : te¥i, t^ + d2^0 (modO}, 

l + cl),^^KMl + l-cl)} = 

llBlif) tfK^Q, 
[Blif] ifK^Q. 

Lfp satisfies condition and if {t, s,p) is a solution to equation arising from 
the newform f then p \ Bi{f). 

Proof. The Proposition follows almost immediately from Lemmas 16.31 IB. 41 □ 

Under the assumptions made (in this Proposition), Method I eliminates all but 
finitely many exponents p provided of course that the integer Bi{f) is non-zero. 
Accordingly, we shall say that Method I is successful if there exists some prime 
I \ 2D so that Bi{f) ^ 0. There are two situations where Method I is guaranteed 
to succeed: 

• If the newform / is not rational. In this case, for infinitely many primes 
I, the Fourier coefficient c; ^ Q and so all the differences ai{Et) — ci and 
I + 1 — ci are certainly non-zero, immediately implying that Bi{f) ^ 0. 

• Suppose that the newform / is rational, and so corresponds to an elliptic 
curve E defined over Q. Suppose that E has no non-trivial 2-torsion. By 
the Chebotarev Density Theorem we know that I + 1 — ai{E) = ftE{¥i) is 
odd for infinitely many primes /. Let I \ 2D be any such prime. ^From the 
models for the Frey curves in Tables ^ 12 El we see that the Frey curve Et 
has non-trivial 2-torsion, and so Z -f 1 — ai{Et) = jJ£'i(F;) is even for any 
value oi t E ¥i, t'^ + d2 ^ 0. In this case ai{Et) — ci — ai{Et) — ai{E) must 
be odd and cannot be zero. Similarly, the Hasse-Weil bound |q| < 2^/1 
implies that / -)- 1 ± q 7^ 0. Thus Bi{f) is non-zero in this case and Method 
I is successful. 



norm map. Lf I \ 2D is prime, let 
B'/if)=lcm{AfK/Q{ai{Et) 



B'Af) 



B'/if) 

lcni{B'/{f),JVK/Q{l 



BiU) 
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8. Eliminating Exponents: Method II 

The second method is adapted from the ideas of Kraus (see also jlSl)- It 
can only be applied to one prime (exponent) p at a time, and if successful it does 
show that there are no solutions to (|10|l for that particular exponent. 

Let us briefly explain the idea of this second method. Suppose / is a newform 
with Fourier expansion as in and suppose p > 7 is a prime. We are interested 
in solutions to equation H1Q|I arising from /. Choose a small integer n so that 
I = np + 1 is prime with I \ D. Suppose {t, s) is a solution to equation (|10|) arising 
from newform /. Then working modulo / we see that d^t^ + D = yP is either or 
an n-th root of unity. Since n is small we can list all such t in Fj, and compute 
ci and ai{Et) for each t in our list. We may then find that for no t in our list are 
the relations in Lemma 16.31 satisfied. In this case we have a contradiction, and we 
deduce that the are no solutions to equation I|1U|) arising from / for our particular 
exponent p. 

Let us now write this formally. Suppose p > 7 is a prime number, and n an 
integer such that I — np + I is also prime and I \ D. Let 

/i„(FO = {CeFr : C" = i}- 

Define 

A{n, /) - |C e ^ln{¥l) : (^^) = or l| . 
For each ( G A{n, I), let 6^ be an integer satisfying 

S'^^ = {C- D)/dl (mod/). 

It is convenient to write a; {Q for a/ (Es^ ) . 

We now give our sufficient condition for the insolubility of (|10|l for the given 
exponent p. 

Proposition 8.1. (Method II) Let D, di, d2 he a triple of integers satisfying (i)- 
(v) of Lemma and p > 7 be a prime satisfying condition 0). Let f be a 
newform with Fourier expansion as in defined over a number field K. Suppose 
there exists an integer n > 2 satisfying the following conditions: 

(a) The integer I — np + 1 is prime, and I \ D. 

(b) Either =-\, or p j AAa7q(4 - c^). 

(c) For all C G A{n, I) we have 

( P\J^K/Q>iai{0-Ci) if 1 = 1 (mod 4), 
\ p\MK/QiaiiCy-cf) ifl = 3 (mod4). 

Then the equation does not have any solutions for the given exponent p arising 
from the newform f . 

Proof. Suppose that the hypotheses of the Proposition are satisfied, and that {t, s) 
is a solution to equation H10() . 

First we show that + c?2 ^ (mod /). Suppose otherwise. Thus + c?2 = 
(mod I) and so / | s. In this case (— f^) = 1, and from (b) we know that p \ 
J^K/qi^ — cf). However, by Lemma we know that ±q = 1 + 1 = 2 (mod *P) for 
some place ^ oi K above p, and we obtain a contradiction showing that t^ + d2 ^ 
(mod I). 
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^From equation H1Q(I and the definition of e in JH)) we see tlie existence of some 
C e v4(n, /) such that 

dlf + D = C, (mod and t = ±5Q (mod 0- 

Replacing t by — < in the Frey curve Et has the effect of twisting the curve by —1 
(this can be easily verified for each Frey curve in Tables^l^lEll- Thus ai{Q — ai{Et) 
if / = 1 (mod 4) and a/(C) — ±ai{Et) if ^ = 3 (mod 4). Moreover, by Lemma ffi. 31 
we know that ai(Et) = ci (mod *P) for some place ^ oi K above p. This clearly 
contradicts (c). Hence there is no solution to (|10|l arising from / for the given 
exponent p. □ 

If the ncwform / is rational and moreover corresponds to an elliptic curve with 
2-torsion, then it is possible to strengthen the conclusion of Proposition 18.11 by 
slightly strengthening the hypotheses. The following variant is far less costly in 
computational terms as we explain below. 

Proposition 8.2. (Method II) Let D, di, d2 be a triple of integers satisfying (i)- 
(v) of Lemma \4-.1\ and p be a prime satisfying condition 0). Let f be a rational 
newform corresponding to elliptic curve E/Q with 2-torsion. Suppose there exists 
an integer n > 2 satisfying the following conditions: 

9 

P 

(a) The integer I = np + 1 is prime, I < — and I \ D . 

(b) Either (^-y^^ = -1; or ai{Ef ^ 4 (mod p). 

(c) For all C G A(n, /) we have 
ai{0^ai{E) if 1 = 1 (mod 4), 



ai{0^±ai{E) ifl = 'S (mod 4). 

Then the equation UfJ\) does not have any solutions for the given exponent p arising 
from the newform f. 

Proof. Comparing this with Proposition 18.11 we see that it is sufficient to show, 
under the additional assumptions, that if a((C)^ = o-i(,EY (mod p) then ai{C,) = 
±ai{E), and if ai{() = ai{E) (mod p) then a;(C) = ai{E). 

Suppose that a;(C)^ = ai{E)'^ (mod p) (the other case is similar). Hence ai(C) = 
±ai{E) (mod p). Now note that both elliptic curves under consideration here have 
2-torsion. Hence we can write 

az(C) = 2bi and ai{E) = 262 

for some integers &i, 62. Moreover, by the Hasse-Weil bound we know that |&i| < 
Vl- Thus 

&i = ±&2 (modp) and I&1 + &2I, |6i - 62I < < p 

p 

since ^ < — . Thus bi = ±62 and this completes the proof. □ 

It remains to explain how this improves our computation. To apply Proposi- 
tion for a particular prime p we need to find a prime I satisfying conditions (a), 
(b), (c). The computationally expensive part is to compute ai{E) = q and ai{(^) for 
all ( € A{n,l). Let us however consider the application of Proposition 18 . 21 rather 
than Proposition The computation proceeds as before by checking conditions 
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(a), (b) first. When we come to condition (c), we note that what we have to check 
is that 

( E^{¥i)^l + l-ai{E) if 1 = 1 (mod 4), 
\ E<;l¥i)^l + l±ailE) if? = 3 (mod 4), 

for each C e A{n,q). Rather than compute ai{Q for each such C, we first pick a 
random point in i?^(F;), and check whether it is annihilated by l + l~ai{E) ifp = 1 
(mod 4) and either of the integers ^ + 1 ± ai{E) if p = 3 (mod 4). Only if this is 
the case do we need to compute ai{Q to test condition (c) in the Proposition. In 
practice, for primes p of about 10^, this brings a 10-fold speed-up in program run 
time for Method II. 

9. Eliminating Exponents: Method III 

Occasionally, Methods I and II fail to establish the non-existence of solutions to 
an equation of the form ()10|l for a particular exponent p even when it does seem 
that this equation has no solutions. The reasons for this failure are not clear to us. 
We shall however give a third method, rather similar in spirit to Kraus' method 
(Method II), but requiring stronger global information furnished by Proposition l3.1l 

Suppose that D, di, d2 are integers satisfying conditions (i)-(v) of Lemma l4. II 
Let Et be one of the Frey curves associated to equation and let / a newform 
of the level predicted by Proposition 16 . 21 with Fourier expansion as in I|1HI . defined 
over a number field K. Define 7/(/) to be the set of r G F/ such that 

• either p \ Nk/q{o-i{Et) ^ ci) and + ^2 ^ (mod I), 

• or p I Nk/q{1 + 1 ± q) and + ^2 = (mod I). 

We suppose that — D is not a square and follow the notation of Section O Fix a 
prime p satisfying ® . Suppose I is a prime satisfying the following conditions: 

(a) l\2D. 

(b) I = np + 1 for some integer n. 

(c) ^ — 7""^ ~ ^' "^^^^ ^ splits in £ = (Q)(^/— £)2) and we let li and [2 be the 

prime ideals above I. 

(d) Each 7 e F is integral at Z; what we mean by this is that each 7 belongs to 
the intersection of the localizations n Oi^. 

We denote the two natural reduction maps by 9i, 02 '■ Cii D Oi^ Fj. These of 
course correspond to the two square- roots for —D2 in F;, and are easy to compute. 

Now let F; be the set of 7 e F satisfying the condition: there exists r G "?!(/) 
such that 

• {diT + DiOii^/^D^Y = 6(1(7)" or (mod /), and 

• {diT + Di92iy/^))" = 6*2(7)" or (mod /). 

Proposition 9.1. (Method III) Letp be a prime satisfying condition @). Let S be a 
set of primes I satisfying the conditions (a)-(d) above. With notation as above, if the 
newform f gives rise to a solution (t, s) to equation Uf)\) . then dit+Di\/ —D2 = 7/3^ 
for some f3 € O and some 7 G n/gsT/. 

In particular, ifCli^s^i is empty, then the newform f does not give rise to any 
solution to equation \1U\) for the given exponent p. 

Proof. Suppose that (i, s) is a solution to equation arising from newform / via 
the Frey curve Et. Clearly 9i{t) = 6*2(0 is simply the reduction of t modulo I. Let 
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T = 6*1 (i) = 92{t) e F;. It follows from Lemma lO that t e Ti{f). Let {x,y) be 
the solution to equation jnj corresponding to {t,s). Thus x — dit. We know by 
Proposition 13.11 that 



for some 7 G F and f3 E O. Applying 6i to both sides and taking n-th powers 
(where we recall that I ^ np + 1) we obtain 



(dir + Did,i^^2)r = e^ilTOm'^^ (mod I). 

However 9i{py~^ = or 1 (mod I). Thus 7 e F; as defined above. The Proposition 
follows. □ 



10. Examples 

It is clear that our three modular methods require computations of newforms of a 
given level. Fortunately the computer algebra suit MAGMA has a package completely 
devoted to such computations; the theory for these computations is explained by 
Cremona |15| for rational newforms, and by Stein 54 5| in the general case. 

Example 2. Absence of Newforms 

Lemma 14.11 and Proposition 16.21 lead us to associate solutions to equation Q 
with p satisfying condition with newforms of certain levels. If there are no 
newforms of the predicted levels, we immediately deduce that there are no solutions 
to equation @. With the help of a MAGMA program we found all values of D in 
the range 1 < < 100 where there are no newforms at the predicted levels. We 
deduce the following result. 

Corollary 10.1. Let D be an integer belonging to the set 

4,16,32,36,64. 

Then the equation ^ does not have any solutions with p satisfying condition 

This Corollary does not add anything new, since equation has already been 
solved by Cohn's method for I? = 4, 16, 32, 36, 64 (but see [HI, g^l, 

Example 3. Corollarv 15.21 solves equation ^ for all values of D in the range ^ 
except for 21 values; these are the 19 values listed in © plus D — 55, 95. As 
indicated in Section [51 the cases D — 55, 95 have been solved by Bennett and 
Skinner. It is however helpful to look at the case D — 95 again as it shows how 
Methods I, III work together in harmony. There is only one possible signature 
(di, c?2) = (1, 95). Thus t = x, s = y and we need to solve the equation 

(12) t^ + 95 = sP, 

under the assumption that p > 7. 

Since di = 1, it follows from Corollarv 15 . 21 that y is even, and so i = a; is odd. 
Replacing t by —t if necessary, we can assume that t = 1 (mod 4). Table pleads 
us to associate the solution {t, s,p) with the Frey curve 
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^From Proposition 16.21 we know that any solution to equation H12|l arises from a 
newform of level 190. Using MAGMA we find that there are, upto Galois conjugacy, 
precisely four newforms at level 190. These are 

h=q-q^-q^+q'-q''+q'^-q' + 0{q% 

f2 = q + - + q" - q'' - iq^ - ^q' + Oiq""), 

f3 = q + q^+q'+q^ + q'+q'-q' + 0(g«), 

^ q - + (t)q^ + + q5 ^ cj>q'^ + <Pq'^ + 0{q^), where + - 4 = 0. 

The first three newforms above are rational, and so correspond to the three isogeny 
classes of elliptic curves of conductor 190. It turns out that none of these elliptic 
curves have non-trivial 2-torsion. By the remarks made after Proposition 17.11 we 
know that Method I will be successful in eliminating all but finitely many exponents 
p. Indeed we find (in the notation of Proposition 17. l|l that B3(/i) = -B3(/3) = 15. 
Thus we know that no solutions to equation (|12(l arise from the newforms /i or /a, 
since otherwise, by Proposition 17. II p \ 15 which contradicts our assumption that 
p > 7. We also find that 53(74) = 2-* x 3 and BriU) = 2"* x 7. Thus no solution 
arises from f^. However, 

S3(/2)=3x7, B7(/2) = 32 X 5 X 7, En(/2) = 0, 

BM) = 3 X 5 X 7 X 13, Bu{f2) = 3^ x 7 x 11. 

We deduce that there are no solutions arising from /2 with exponent p > 7. It does 
however seem likely that there is a solution with p — 7. Moreover, an attempt to 
prove that there is no solution with p — 7 using Method II fails: we did not find 
any integer 2 < n < 100 satisfying the conditions of 18. II 

We apply Method III (and follow the notation of SectionEJ. Write uj = 
Taking 

5" = {113,127,239,337,491} 

we find that 

r-528- 2cj' 
""'^'^^ - 1 2187 
Thus if we have any solutions at all then, by Proposition 19. II we know 

for some integers U, V. Equating imaginary parts and clearing the denominators 
we find that 

- 2834013V"^f/^ - 29U37W^U + 2818247^^ = 2187. 

Using pari/gp we find that the only solution to this Thue equation is given by 
[/ = -3, y = 0. This shows that {t, s) = (529, 6). 

The reader will notice that {t, s) = (—529, 6) is also a solution to equation l|12|l 
with p ^ 7 but it seems to have been 'missed' by the method. This is not the 
case; we are assuming that the sign of t has been chosen so that t = 1 (mod 4). 
The solution (t, s) — (—529, 6) arises from some other newform (probably at some 
different level) and via a different Prey curve which we have not determined. 
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Example 4. For our last example we look at the case where D — 25. This, like 
18 other cases, must be resolved by a combination of the modular approach and 
our lower bound for linear forms in three logarithms which is to come. We assume 
that p > 7, and so p satisfies conditions ((BJ. There are now two possible signatures 
((ii,c?2) — (1,25), (5,1) satisfying the conditions of Lemma l4.ll However, by 
CoroUarv 15.21 we may suppose that di > 1 and so di = 5, (i2 = 1. We write 
t = a;/5, s = y/b where we know that s are integral by Lemma EITI Equation 
becomes 

+ 1 = 5P-2sf , t ^ ±1. 

Following Tabled we associate with any solution to this equation the Frey curve 

Et: Y"^ =X'^ + 2tX^ - X, 

and we know by Proposition 16.21 that any solution must arise from a newform of 
level 160. Using the computer algebra system MAGMA we find that there are upto 
Galois conjugacy three such newforms: 

h=q-2q^-q'-2q' + 0{q''), 

f2=q + 2q^-q^ + 2q' + 0{q^), 

h = q + 2V2g3 + g5 _ 2V2g^ + 0{q^). 

The first two newforms are rational, corresponding respectively to elliptic curves 
160A1 and 160B1 in Cremona's tables il5j. The third has coefficients in i^T = (Q(\/2) 
and is straightforward to eliminate using Method I. In the notation of Proposi- 
tion l7.1l we find that if /a does give rise to any solutions (i, s,p) thenp | i?3(/3) = 24. 
This is impossible as p > 7, and so /s does not give rise to any solutions. 

We where unable to eliminate newforms /i , /2 using Method L Instead using our 
implementation of Method II in MAGMA we showed that there are no solutions arising 
from either form with 7 < p < 100. With our implementation of the improved 
Method II (Proposition IS . 2|l in pari/gp we showed that there are no solutions with 
100 < p < 163762845; this took roughly 26 hours on 2.4 GHz Pentium IV PC. The 
choice of where to stop the computation is of course not arbitrary, but comes out of 
our bound for the linear form in logarithms. We will later prove that p < 163762845 
thereby completing the resolution of this case. 

11. Results III 

We applied the methods of the previous sections to solve all equations Q with 
D is our range (0). 

Lemma 11.1. Suppose D is in the range (0) and p is a prime satisfying 
Suppose (x,y,p) is a solution to equation 0) that is not included in the tables. 
Then D is one of 

(13) 7, 15, 18, 23, 25, 28, 31, 39, 45, 47, 60, 63, 71, 72, 79, 87, 92, 99, 100. 

Moreover {x,y,p) has signature (^1,^2) and arises from an elliptic curve E and 
p > Po where E, po and (^1,^2) are given by Table^ 

Proof. We wrote a MAGMA program that does the following: For each D in the 
range (O we write down the set of possible signatures (c?i, o?2) satisfying the condi- 
tions of Lemma 14.11 
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Table 4. Computational details for Lemma [11.11 and its proof. 



D 


(^1,^2) 


E " 


Po 


Machine " 


Time 


7 


(1,7) 


14A1 


181000 000 


PI 


26h, 43mn 


15 


(1,15) 


30A1 


624 271 465 


SI 


252I1, 50mn 


18 


(3,2) 


384D1, 384A1, 
384G1, 384H1 


306111726 


S3 


293I1, 14mn 


23 


(1,23) 


46A1 


855 632 066 


S2 


477I1, 36mn 


25 


(5,1) 


160A1, 160B1 


163 762 845 


P2 


25h, 58mn 


28 


(2,7) 


14A1 


315 277186 


PI 


55h, 41mn 


31 


(1,31) 


62A1 


860111230 


S3 


2421i, 2mn 


39 


(1,39) 


78A1 


852 830 725 


PI 


I93I1, 41mn 


45 


(3,5) 


480B1, 480F1, 
480G1, 480H1 


340 749 424 


SI 


4481i, 43mn 


47 


(1,47) 


94A1 


1555 437629 


S3 


45II1, 34mn 


60 


(2,15) 


30A1 


358 541296 


SI 


I3OI1, 30mn 


63 


(1,63) 


42A1 


292 825 735 


SI 


99h, 45mn 


71 


(1,71) 


142C1 


2 343 468 548 


S3 


6971i, 26mn 


72 


(3,8) 


96A1, 96B1 


451620 034 


SI 


3I6I1, 27mn 


79 


(1,79) 


158E1 


1544 381661 


S3 


4481i, 47mn 


87 


(1,87) 


174D1 


1 148 842108 


S3 


329I1, 45mn 


92 


(2,23) 


46A1 


996 255151 


S3 


2851i, lOmn 


99 


(3,11) 


1056B1, 1056F1 


593 734 622 


P2 


1381i, 46mn 


100 


(5,4) 


20A1 


163 762 845 


PI 


21h, 23mn 



"We give here the Cremona code for the elliptic curves E as in his book il5i and his online tables; 
Kittp : //www . maths . nott . ac . uk/p ersonal/jec/f tp/data/INDEX-html] 
The machines are as follows 

PI 2.2 GHz Intel Pentium PC. 

P2 2.4 GHz Intel Pentium PC. 

51 Dual processor 750MHz UltraSPARC III. 

52 650 MHz UltraSPARC He. 

53 UltraSPARCIII with 12 processors of 1050 MHz speed. 



For each such pair (di, ^2) write down the (one or two) Frey curves given by the 
Tables n] 12 O bearing in mind the information given by Corollarv l5.2l 

For each Frey curve we compute the conductor (given by Proposition 16. 2|) of the 
newforms giving rise to possible solutions, and then write down all these newforms. 

We attempt to eliminate each newform / using Method I. This involves searching 
for primes / \ 2D such that (in the notation of Proposition 17. l|l Bi{f) ^ 0. If we 
are successful and find such primes /i, . . . , then by Proposition l7. li the exponent 
p divides all the Bi. (/), and so divides their greatest common divisor B (say). If B 
is divisible by any prime p that satisfies condition © then we attempt to eliminate 
this possible exponent p using Method II: this involves searching for an integer 
2 < n < 100 satisfying conditions (a), (b), (c) of Proposition 18.11 If one such n 
is found then we know that there are no solutions for the particular exponent p. 
Otherwise we apply Method III (Proposition I9.1|l to write down Thue equations 
leading to possible solutions (see below). 
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As predicted by the comments made after Proposition 17. II Method I succeeded 
with aU non-rational newforms and all rational newforms corresponding to elliptic 
curves with only trivial 2-torsion (it also succeeded with some rational newforms 
corresponding to elliptic curves with non-trivial 2-torsion). Indeed, we found no 
solutions arising from non-rational newforms for D in our range |2| 

We are left only with rational newforms / that correspond to elliptic curves E 
having some non-trivial 2-torsion. The details of these are documented in Table 01 
For primes p < 100 satisfying condition © we attempt to show that there are no 
solutions arising from E for the particular exponent p using Method II (as before). 
If this fails for a particular exponent p, then we use Method III to write down the 
Thue equations leading to the possible solutions. 

Our proof that p > 100 is now complete except that there are some Thue equa- 
tions to solve. We had to solve Thue equations of degree 7 for D = 7, 47, 79 and 95. 
These were solved using pari/gp and the solutions incorporated in our Tables. We 
also had to solve a Thue equation of degree 11 for D = 23, of degree 17 for D = 28, 
and of degree 13 for D = 92. We were unable to (unconditionally) solve these three 
Thue equations using the in-built functions of pari/gp. The reason is that, in each 
case, it was impossible for pari/gp to prove that the system of units it had found 
- though of correct rank - was maximal. We are grateful to Dr. Guillaume Hanrot 
for sending us his pari program for solving Thue equations without the full unit 
group. This program, based on '19', solved all three equations in a few minutes. 

For the next step we implemented our improved Method II (Proposition I8.2|l 
in pari/gp (see the remark after the proof). To complete the task and show 
that p > pq for any missing solution we used our pari/gp program to disprove the 
existence of any missing solution for each prime 100 < p < po- We ran this pari/gp 
program on various machines as indicated in Tabled The total computer time for 
this step is roughly 206 days. □ 

Remark. The reader may be surprized that some of the computations were done 
in MAGMA while others were carried out in pari/gp. As stated earlier, MAGMA has a 
package for computing modular forms. This is essential for us, and is unavailable 
in pari/gp. 

For showing that p > po, it is simply not practical to use MAGMA. Here we are 
using the improved Method II ('Proposition 18 . 2|l . The main bottle-neck in Method 
II is computing ai(E) for primes I that can be about 10^^ (recall I is a prime 
satisfying I = 1 (mod p)). For this pari/gp uses the theoretically slower Shanks- 
Mestre method ^21 rather than the theoretically faster Schoof-Elkies-Atkin j^HI 
method used by MAGMA. But for primes of the indicated size it seems that pari/gp 
is about 10 times faster than MAGMA. 

The reader may also note that two of the machines we used are multiprocessor 
machines. The computation for each D could have been speeded up considerably by 
parallelising. We however decided against this, so as to keep our programs simple 
and transparent. 

12. The 'Modular' Lower Bound for y 

In this section we would like to use the modular approach to prove a lower 
bound for y with D in the range (jSJ. Before doing this we prove a general result 
for arbitrary non-zero D. 
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Proposition 12.1. Suppose D is a non-zero integer, and di, d2 satisfy (i)-(v) of 
Lemma \4-1\ Suppose {t, s,p) is a solution to equation HfJ\) arising from a rational 
newform f via a Frey curve Et- Then either rad(s) | 2di or \s\ > {^/p — 1)^. 

Proof. Since the newform is rational we know that the newform / corresponds to 
an elliptic curve E/Q whose conductor equals the level of /. 

Suppose rad(s) does not divide 2di. Since t and d2 are coprime we see that there 
is some prime Z | s so that I \ 2D. By Lemma 15^ we see that p divides l + l±ai{E). 
It follows from the Hasse-Weil bound that I + 1± ai{E) ^ 0, and so 

p< l + \±ai{E) < {Vl + lf, 

where again we have used Hasse-Weil. Thus I > (y/p — 1)^. The Proposition follows 
as Z I s. □ 

Corollary 12.2. Suppose D is one of the values in If{x,y,p) is a solution 

to equation ^ not in the Tables below then y > (y^ — 1)^. 

Proof. Suppose D is in the range ijSJl and (x,y,p) is some solution to equation @ 
not in the below Tables. From the preceding sections we know that this solution 
must satisfy condition JB)). Moreover by Lemma l4.1l 

X — dit, y — rad(di)s, 

where {t,s,p) satisfy equation (|10|) for some di, c?2 satisfying conditions (i)-(v) of 
that Lemma. 

We have determined for D in the specified range all solutions to equations Hl()|l 
arising from non-rational newforms (indeed there were none) . Thus we may suppose 
that our putative solution arises from a rational newform. By Proposition 1 1 2 . l1 we 
see that either \y\ > \s\ > {^/p~■ 1)^ or rad(s) | 2di. We must prove that the second 
possibility does not arise. 

Suppose that rad(s) | 2di. From Lemma EH we see that rad(y) | 2di. We first 
show that rad(y) ^ 2. For in this case we have reduced to an equation of the form 

+ D = 2"^. For |-D| < 2^^, which is certainly the case in our situation, Beukers 
01 Corollary 2] shows that 

m< 18 + 21og|D|/log2. 

A short MAGMA program leads us to all the solutions to this equation for 2 < D < 100 
and we find that these are already in our tables. 

Thus we may suppose that rad(y) | 2di and rad(y) ^ 2. An examination of the 
possible cases reveals the following possibilities 

(1) D = 18, 45, 72, 99 and rad(?/) = 3, 

(2) D = 25, 100 and rad(y) = 5. 

On removing the common factors, each case quickly reduces to an equation that 
has already been solved. For example, we must solve x'^ + 100 = under the 
assumption that rad(t/) = 5 or equivalently the equation x'^ + 100 = 5'". Removing 
the common factor reduces to the equation + 4 — S™^'^. But the equation 
+ 4 = y" has already been solved and has only the solutions {X, Y, n) ~ (2, 2, 3), 
(11,5,3). We quickly see that the only solution to x^ + 100 = y^ with y = 5 is 
(2^, 2/, P) = (55,5,5). □ 
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13. The Linear Form in Logarithms 



It is useful at this point to recap what we have done so far. We would like to 
complete the proof of Theorem ^ by showing that our Tables at the end are not 
missing any solutions. So let us suppose that our Tables at the end are missing 
some solution {x,y,p) to equation Q for some value of D in our range (01. We 
have proved fLemma lll.l(l that D is one of the values in p3|l . Moreover, (again by 
Lemma 111 . II and by Corollarv ll2.2|l any missing solution (x,y,p) must satisfy 

(14) P>Po, y>(VP-i)'> 

with po given by Table 01 Our aim is to show that p < po thus obtaining a 
contradiction. 

^From the table of values of po we know that 

(15) y, p>108 

and indeed much more, though this inequality is sufficient for much of our later 
work. In the remainder of this paper we assume that D is one of the remaining 
values (|13|l . and always write (as before) D = D\D2, where D2 is square- free. The 
triple {x,y,p) will always be a solution to equation Q supposedly missing from 
our Tables and hence satisfying the above inequalities. 

In this section we write down the linear form in logarithms corresponding to 
the equation Q and apply a Theorem of Matveev to obtain upper bounds for the 
exponent p. These upper bounds obtained from Matveev's Theorem are not small 
enough to contradict our lower bounds for p obtained in Lemma 111. II but they are 
needed when we come to apply our bounds for linear forms in three logarithms 
given in the next section. 

Lemma 13.1. Let (^1,^2) be the signature of our supposedly missing solution 
{x,y,p) (which we know from Lemma \TJ~l\l . Define 

^^Ui, forD^l (mods), 
\2di, forD = 7 (mods). 

Then d is a prime power, say d — q'^ for prime q, where moreover, q splits in 
C = Q(\/~-D2), say (q) — qq. Let ko be the smallest positive integer such that the 
ideal q'^" is principal, say q*^" = (ao). Also let 

\ko, if ko is odd, J 2, if ko is odd, nko 

k = < and k, ~ < so that k — ——. 

\ko/2, if ko is even, II, if ko is even, 2 

Then there exists 7 G £ such that 



x-Diy/^D2 
x + Di^/^D^ 



where 

a = ao/ao, h{a) = , h(7) = — - — . 

K I 

Proof. We begin with the factorization 



ix + Diyf^2Kx-Di^^2)^y''- 

Our first step is to show that any prime divisor q oiy splits in C. Suppose otherwise, 
then we may write {q) = q or (q) — q^ for some prime ideal q satisfying q = q. If 
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p — 2r + 1 then clearly q*" divides both factors on the left-hand side above, and so 
divides 2Di\/ —D2. This is impossible in view of the fact that p is enormous, and 
2 < Z? < 100. Thus we have shown that every prime divisor q of y splits in C. 
Let us write 

y^^JIft"' and (<7i) = q, 7^ q^, i G /. 

Then 

iei 

where we assume (for commodity of notation) that bi > Ci for all i, and thus 



iei 



with 
Let 

clearly 



bi + Ci = pai, for all i £ I. 
= gcd (x + Di^-D2,x- I?i 



This shows that O = {d) where d G Z. We would like to calculate this d and verify 
that its value is in agreement with pt)|l . /^From the definition of we see that d \ 2x 
and d \ 2Di. However, by our definition of signature, gcd(x^,£') c?^. It follows 
that d^ I Adi and so d \ 2di. But di \ x and di \ Di. Hence di \ 5 and so di \ d. 
Thus d = di or d = 2di. 

We note the following cases: 

• U D2^7 (mod 8) then 2\y. Thus 2 f d and so d = di. 

• Suppose Z?2 = 7 (mod 8) . Now from Lemma 14.11 and its proof we know 
that D = d\d2 and x — dit where gcd(i, d2) = gcd(di,d2) — 1. Clearly 
d2 = d§£'2 with da = Di/di integral. Suppose first that di is even. It 
follows easily that t, d2 are odd and 



(d)^o^2d,r+"^^,^-"^^ 



Hence (2di) | d and so d = 2di. 

The only case left to consider is D2 = 7 (mod 8) and di is odd. By exam- 
ining Tabled we see that di = 1. Thus 2 | y by Corollarv l5.2l Clearly x is 
odd, and the same argument as above shows that d = 2 = 2di. 
This proves that d satisfies (llfill . By looking again at the possible values of di in 
Table 21 we see that d is a prime-power in all cases. Let j E I such that d ~ q^^ 
Thus c, = for all j ^ i. Then 

(x + i?iv^) = q;^-q^-ni''"'' 

whence 

(x + D,^^2) - (q, qr')'^ ■ n l^"""' = ^"') 0"' 
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where a and q are integral ideals with 

o = q;^ AA(a) = g/^ = rf, Mig) ^ y, 
and jV denotes the norm. Thus, as ideals, 



We define feg, k, k, ao as in the statement of the Lemma. Thus a'^" = (ag) and 
we have the relation (between ideals) 



However p is an enormous prime certainly not dividing the class number. This 
shows that is also principal, = (70), say, where 70 is an algebraic integer 
chosen so that the following equality of elements of C holds 



Note that 

AA(ao)=rf''°, AA(7o) = /. 

Write 

a = ao/ao, 7 = ±70/70. 

The proof of the Lemma is complete except for the statements about the heights 
of a, 7. These follow from Lemma [13.21 below. □ 

Lemma 13.2. Let a he an algebraic number whose conjugates are all (including a 
itself) 0/ modulus equal to 1, then 

= logo, 
deg a 

where a is the leading coefficient of the minimal polynomial of a. In particular, if 
a — ao/o^o where ao is a non-real quadratic irrationality, then 

h(a) = ilogAA(ao). 

Proof. Let d = deg a. By hypothesis a is a root of a polynomial of 1\X] of the 
form P{X) = aX'^ + ■■■. We have 

h(a) = ilogM(P), 

where M denotes Mahler's measure, and the first result easily follows since the roots 
of P are of modulus equal to 1. This proves the first assertion. 

The proof of the second assertion, which is quite easy, is omitted. □ 



We now write the linear form in three logarithms. Define 

A = logf-"^^^ 

where we have taken the principal determination of the logarithm. 
Lemma 13.3. 

log|A| < -| log2/ + log(2.2Di/D^). 
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Proof. We will rely on the lower bounds (|15|l . Clearly 



< 2- 



A standard inequality (01], Lemma B.2) shows that 



|A| < 2.1 ■ 

so that 



log|A| <-\og\x\+\og{2.lDi^\D2\). 

Using the fact that — — D, and a similar argument to the one above, we 
deduce the Lemma. □ 

The main tool to bound p will be the theory of linear forms of (at most three) 
logarithms. We need the special case of three logarithms of the Theorem of E. M. 
Matveev. 

Theorem 2 (Matveev). Let Xi, A2, A3 be Q-linearly independent logarithms of 
non-zero algebraic numbers and let bi, 62, 63 be rational integers with 61 7^ 0. 
Define aj = exp(Aj) for j = 1, 2, 3 and 

A ^ biXi + 62A2 + b^X^. 

Let V be the degree of the number field Q{ai, a-i, 0^3) over Q. Put 

X= [M(ai,a2,a3) : K]. 

Let A\, A2, A'i be positive real numbers, which satisfy 

Aj > max{Vh{aj), |Aj|,0.16} (1 < j < 3). 

Assume that 



Define also 



Then 



B > max|l,max{|6j|Aj/Ai; 1 < j < 3}|. 
= ^^e^(7 + 2x) (|)'(20.2 + log(3^-^P^log(eP) 



log I A| > -Ci Ai A2A3 log (1.5 eVB log(eX')) . 
For V — 2 and x — 2, this gives 

(17) log |A| > -1.80741 X 10"yli A2A3 log (13.80736 B), 
whereas, for V — 2 and x = 1, we get 

(18) log I A| > -7.25354 x lO^°Ai A2A3 log (13.80736 B). 

Proof See □ 
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13.1. A Preliminary Bound for p. It follows from Lemma Fl 3. II that 

kA = Kloga +plog7 + im = reloga +plog7 + r log(— 1), 

which appears as a linear form of logarithms. But a small transformation of this 
form leads to better estimates. Write 

kA — Klog(£ia) +plog(e27) + 

where Si and e2 = ±1 are chosen in such a way that 

|log(eia)| < tt/2 and |log(£27)l < 7r/2, 

where we take the principal values for the logarithms, and where r e Z is such that 
|A| is minimal (we keep the same notation r as before for simplicity). 

Remark. Indeed, we can take any roots of unity in C for ei and £2- The only relevant 
case for our set of outstanding values of D are D = 25, 100, where C = 
whence we can realize 

I log(£ia)| < 7r/4 and | log(£27)| < 7r/4, 

and we write 

A = 2 log a+p log 7 + r log (, 

where ( = e*'^/^. 

We now return to the general case. By Lemma [13. 31 

log|fcA| < -| logy + \og{2.2 kDi^). 

Our lower bound for x, y, p implies that log|/cA| is very small and it is straight- 
forward to deduce that 

I I ^ 

I I - 2 

We can write fcA in the form 

kA = biXi + 62A2 + 

with 

61 = K (= 1 or 2), ai — £iq;, 62 — p, a2 — £27, ^3 — ~ -1 

and 

= -logd, Ai=logai, h(a2) = — - — , IA2I < 7r/2, h(a3) = 0, A3 = ztt, 

K I 

except for the case L — 1) studied in the previous remark where A3 = m 12. 

Applying Theorem 12 we have T) — x = ^ ^-nd we can take 

^ f2fclogd ttI ^ r, , TTT 

Ai=:max<^ ' r ' ^2 = max|fclog?/, -| , A3 = tt 

and (using some change of numcrotation in Theorem 12) 

S =p+ 1 

(this choice of B is justified by the inequality \r\ < (p + l)/2 proved above), and 
we get 

p<C2k^ log(2i?i)logp. 

This implies 

V < Csk"" \og{2Di) log{k^ log(2Z?i)) , 
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and thus 

(19) p < CiD2 log(27^i) log(i?2 log(2i?i)) , 

where the constants could easily be explicated. 

Lemma 13.4. Suppose D is one of the remaining values H,^) and {x,y,p) is a 
solution to 0) missing from our Tables. 



• If D^l thenp< 6.81 x lO^^. 

• Otherwise if D is square-free then p < 1.448 x 10^^. 

• For other values of D, we have p < 3.966 x 10^^. 

15 



Thus in all cases p < 1.5 x 10 

Proof. This is a simple application of Matveev's Theorem [3 If D = 7 it is easy 
to show that the ao arising in Lcmma ri3.il is (upto conjugation) (1 + \/—7)/2, we 
know that fc = 1; thus 7V(ao) = 2 and ^(logao) = 1.2094292028 . . . Then we can 
apply Theorem 121 with £> = 2, x = 2 and 

Ai=7r/2, = logy, logAg^TT, B=p+1. 

After a few iterates we get the stated bound on p. 

In the application of Theorem|3 we can take, for all the squarefree values of D, 

J 7 log 2, if fco is odd, J 7 logy, if ko is odd, 

I 8 log 2, if ko is even, I 4 logy, if ko is even, 

so that 

A1A2 < 49 log 2 X logy 

and we get 

p < 1.448 X 10^^ 
For all the remaining values of Z?, we can take 

'logic, 




iih = 


1, 




flogy. 


iih = 


1, 




if /l = 


2, 




1 logy. 


iih = 


2, 


A3 = TT 


if = 


3, 




[3 logy. 


iih = 


3, 





A1A2 < 9 log 2 X logy 
p < 3.966 X IQi'*. 

14. A NEW ESTIMATE ON LINEAR FORMS IN THREE LOGARITHMS 



□ 



We present the type of linear forms in three logarithms that we shall study. 
We consider three non-zero algebraic numbers ai, a2 and 03 and positive rational 
integers 61, 62, ^3 with gcd(&i, 62, ^3) = 1, and the linear form 

A = 62 log Q!2 - h log ai - bi log ^0. 

We restrict our study to the following cases: 

• the real case: ai, a2 and are real numbers > 1, and the logarithms 
of the ai are all real (and > 0), 
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• the complex case: ai, a2 and as are complex numbers of modulus one, 
and the logarithms of the ai are arbitrary determinations of the logarithm 
(then any of these determinations is purely imaginary). 

In practice this restriction does not cause any inconvenience since 

|A|>max{|3?(A)|,|3(A)|}, 

and so we can always reduce to the above cases. 

Following [2], we use Laurent's method, and consider a suitable interpolation 
determinant A. 

Without loss of generality, we may assume that 

62I loga2| = 6i| logail + 63] log 03] ± |A|. 

We shall choose rational positive integers K, L, R, S, T, with K, L > 2, we 
put TV = K'^L and we assume RST > N. Let i be an index such that {ki,mi,ii) 
runs trough all triples of integers with Q < ki < K — 1, Q < rrii < K — 1 and 
Q < ii < L — 1. So each number 0, . . . , K — 1 occurs KL times as a hi, and 
similarly as an m^, and each number 0, . . . , L — 1 occurs times as an £i. 

With the above definitions, let 



A = det 



II I CIt CI9 Cl-^ 



where rj, Sj, tj are non-negative integers less than R, S, T, respectively, such that 
(rj, Sj,tj) runs over N distinct triples. 
Put I3i =bi/b2, fiz^bi/b2. Let 

n L-1 R-1 ^ S~l ^ T-1 ^ S~l 

and 

4 

(K-l \ K{K~l) 
n ^-j 

Following 15^1, Lemme 8, we can prove that 

^ (i?-l)b2 + (g-l)&l , , {T-l)b2 + {S~l)bs 

log b < log ^ — ^ + log 

21og(2^j^/e3/2) 2 + 671-2 + log g 

— Z log A + O ; . 

^ K-l 3K{K-l) 

Then, we have X^iLo^ '^^ ^ ^ ^^"^ (.2.' formula (2.1)) 

fir, £iSj fitj _ X^irj+Sj/li) \,{tj+Sjf33)f^ o ./n 

"1 "2 "3 - Q!i as [L + OtjA), 

where 



A' = |A| • max 
and where all \9ij\ are < 1. 



fLiJe^^^l^l/t^&i) xS'e^-5|'^l/(2fc2) ^2^gLT|A|/(2fc3) 
1 2\h\ ' 2\h\ ' 2N 
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14.1. An upper bound for |A|. Put 

N N N 



i=i j=i i=i 



and 



_ 1 N _ NLR _ NLS ^ _ NLT 

^-I'URSf' ^'-^^^ ^'-^^^ ^'-^3, 

then, see 0: 

Proposition 14.1. With the previous notation, if K >3, L > 5 and A' < p 

for some real number p > 1, then 

log I A I < ^ M, log I a, I + p ^ G, I log a, I + log(iV! ) + log 2 + - (if - 1 ) log 6 
/NKL NKL NK N\ , 

14.2. A lower bound for |A|. Using a Liouville estimate as in [2^1 Lemme 6, we 
get (see 0): 

Proposition 14.2. If A ^ then 

log I A I > - N\ogN + J2 + ) log I a, | 

3 

- 2D ^ G,h(aO - ^^-^ {K-\)N log b. 

i=l 

14.3. Synthesis. We get (see again 

Proposition 14.3. With, the previous notation, if K > 3, L > 5, p > 1, and if 

A 7^ then 

A' > p-^^ 

provided that 

+ D(A:- 1) log6-21og(e/2), 

where the ai are positive real numbers which satisfy 

ai> p\ log ai\- log \ai\+2Dh{ai), i = 1, 2, 3. 

Remark. We notice that the statement of Proposition ll4.3| is perfectly symmetric 
with respect to the 6i's or the a^'s, except for the choice of b. From now on we do 
not assume that bi and 63 are positive, but we still suppose that 62 > and that 

62I loga2| = \bi logail + I63 log as] ± |A|. 
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14.4. A zero-lemma. To conclude we need to find conditions under which one of 
our determinants A is non-zero, a so-called zero-lemma. We use a zero-lemma due 
to M. Laurent ii24i which improves "18' and provides an important improvement on 
the zero-lemma used in our previous paper 9 : 

Proposition 14.4 (M. Laurent). Suppose that K , L are positive integers and that 
Si, Yi2 and S3 are finite subsets 0/ x C* containing the origin and such that 



(i) 
and 

(ii) 



Caid{\xi+ ^1X2 : {xi,X2,y) £^1} > K, V(A, ^) 7^ (0, 0), 
Card{2/ : {xi,X2,y) G Si} > L, 



(M) 



Card{(Axi + ^X2,y) : ixi,x2,y) € S2} > 2KL, V(A, fi) ^ (0, 0), 
Card{(xi,a;2) : (a;i, a;2, y) G S2} >2K^, 
and also that 

(iii) Card S3 > eiCL^ 

Then, the only polynomial P G C[Xi, X2,Y] with deg_^. P < K for i — I, 2, and 
degy P < L which is zero on the set Si + S2 + S3, is the zero polynomial. 

We now study the above conditions in detail. For j — 1, 2, 3, we shall consider 
finite sets Sj defined by 

Sj- = {(r + s^i,t-f s^3,a>^4) : < r < Rj , < s < Sj , < t < T^} 

where Rj, Sj and Tj are positive integers and where 

R -^^ R -^^ 

Pi — 7-, P3 — T-- 

02 02 

In practical examples, generally the following condition holds: 

{either ai, a2 and are multiplicatively independent, or 
two multiplicatively independent, the third a root of unity 7^ 1. 

We also assume that 

(11) Card{(a:i,X2) : (xi, X2, zj) G Si } = (i?i + l)(5i + l)(ri + 1), 
and 

(12) Card{(xi,X2) : (xi, 2:2, y) G S2} = (i?2 + 1)(^2 + l)(r2 + 1). 
Notice that if 

(r -I- s/3i, t -I- s/?3, a^ajag) = (/ + s'/3i, t' + s'/?3, Q^i ^2 "3 ) 

then, when hypothesis (M) holds, two pairs of the integers (r, s, t) and (r', s' , t') are 
equal which clearly implies that indeed these triples are equal: for example if ai and 
a2 are multiplicatively independent, then the equality a^^ajOfg — Oi\ implies 

r = r' and s — s' and then we conclude that t = t' (use the relation X2 — X2). 
Hence 

(M) => CardSj- = {R, + l){Sj l)(Tj + 1), j = 1, 2, 3. 

The conditions of the zero-lemma are the following: 
(i) The first condition is divided into two subconditions 

(i.l) Card{Aa;i +^a;2 : (a;i,a;2,y) G Si} > is:, V(A, ^i) 7^ (0, 0). 
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This is the most technical of the above conditions, we study it in detail later. 
The second subcondition is 

(i.2) Card{?/ : {xi,X2,y) e Si} > L. 

Wc use now hypothesis (M), and we also notice that the order in C* of a root of 
unity ^ 1 is at least equal to 2 (since ^ 1), thus this condition is satisfied if 

2(i?i + + l)(Ti + 1) 



(C.i.2) 
where Wi is defined by 



>L, 



' Ri, if ai is a root of unity. 

Si, if a2 is a root of unity, 

Ti, if as is a root of unity, 

1, otherwise. 



But see also the remark after (C.ii.l) below. 

(ii) The second condition of the zero-lemma is also divided into two subconditions, 

the first being 

(ii.l) CsiTd{{\xi+ ^lX2,y) : {xi,x2,y) € S2} > 2KL, V(A,m) ^ (0,0). 

We replace it by the stronger condition 

Cardjj/ : (.xi,.T2,?/) G S2} > 2KL. 

Then, by the study of the case (i.2), we see that it is enough to suppose that (when 
condition (M) holds) 

(i?2 + l)(S2 + l)(T2 + l) 



(C.ii.l) 

where W2 is defined by 



W2 + I 



> KL, 



W2 



i?2, if a 1 is a root of unity, 

^2, if Q!2 is a root of unity, 

T2, if 0:3 is a root of unity, 

1, otherwise. 



Remark. When (for example) as is a root of unity of order v, condition (C.ii.l) can 

be replaced by 

(C'.ii.l) ly (i?2 + 1){S2 + 1) > 2KL, 

and condition (C.i.2) can be replaced by 

(C'.i.2) u {Ri + l)iSi + 1) > L. 

The second subcondition of condition (ii) of the zero-lcmma is 
(ii.2) Card{(xi,a;2) : {xi,X2,y) G S2} > 2K'\ 

By (I2) this condition is equivalent to 
(C.ii.2) (i?2 + 1)(S2 + 1)(T2 + 1) > 2K^. 

(iii) There is just one condition, namely 

Card S3 > 6KL'^. 
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When (M) holds, this is equivalent to 

(C.iii) {R3 + l){S3 + l){n + l)>6K^L. 

Now we have 'translated ' all the conditions of Proposition 114.41 except the 
subcondition (i.l). We come back to this situation in the following Lemma which 
brings some extra information to Proposition 3.1.1 of 

Lemma 14.5. Let A, B and C be non-zero rational integers with gcd(yl, B,C) = 1 
and let D he an integer. Define 

n = {{x,y,z) e : Ax + By + Cz^D] 

and consider the set 

E = {{x,y,z) <El? : Q <x < X, Q <y <Y, Q < z < Z], 

where X , Y and Z are positive integers. Let 

M = Card{(x,?/,z) G S : Ax + By + Cz^ D]. 

Then 



M < I 1 

where 



1 



Y 



\C\la 



id M < [I 





X 




Z 








)(- 




) 




a 




[\B\/a\ 





a = gcd(B,C). 

// we suppose that 

M > max{X + Y + 1,Y + Z + 1, Z + X + l} 

then 



\A\< 
\C\< 



(y + i)(z + i) 

M - max{r, Z} ' 
(X + l)(y + l) 
M - max{X, Y} ' 



\B\ < 



(X + 1)(Z+1) 
M - max{X, Z} ' 



Proof. If the image (by the map (a;, y, z) i~> Ax + By + Cz) of a point (x, y, z) G Z"^ 
belongs to the plane 11 then 

Ax = D (mod a), 

where A and a are coprime since gcd{A, B, C) = 1. This shows that the number of 
such X which satisfy < x < X is 

X 

a 



< 1 



To simplify the notation we suppose for a while that A, B and C are positive. Let 
now x be fixed, with < x < X , and such that the images of two points {x, y, z) 
and {x,y',z') belong to 11. Then 

B{y'-y) = C{z^z'), 

where we suppose (as we may) that y is minimal (then y' > y). Hence there exists 
fc G N such that 

y'~y — k(C/a) and z — z'~k{B/a). 
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It follows that, for x fixed, the number of (x, y, z) G S whose image belong to 11 is 

Y 



Hence 



M < 1 



< 1 + 



C/a. 



1 



Y 



which proves the first inequality of the Lemma. The proof of the second one is the 
same (looking at z). 

For ^ > 1 put 



X 



then 

Suppose now 
Put 



M</(q). 
M > max{X + 1,Y +1, Z 



ai = max{l, C/Y}, a2 = min{C, X}. 

• If C > y and 1 < a < C/Y then we get M < X + 1, contradiction, thus 

(XY 
1 + — 

• If C > X and a > X then wc get M < Y + 1, contradiction, thus 

XY\ 



OX 



a < a2 and f{a2) =2(1 + 



If C < min{X, Y} then ai = 1 and a2 = C and 

Y 



C 



/(ai) = (X + l) 1 + - , f{a2)= 1 + - (F + l). 



X 



C 



It is easy to check that /" is positive and, from the previous study, it follows 
that 

M<max{/(ai),/(a2)}. 

Considering the different cases C > max{X, K}, X<C<Y,Y<C<X and 
C < min{X, Y} we get always 



Y 



M < max<^ (X + 1) 1 + - , 1 + - (y + 1) 



C 



X 



UX >Y then 



C 



M < (X + 1) 1 



'(x + i)(i + ^), ifx>y, 

(1 + (y + 1), otherwise. 



y 

Cj ' 



which implies 



M-iX + l)<^^^, hence C< + 



C 



M- (X + 1)' 
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and the hypothesis M > X + Y + I leads to 

(x + i)(y + i) 

otherwise {i.e., if X < Y) we get 

(x + i)(y + i) 

Finally, we always have 

{X + 1){Y + 1) 



\C\ < 



M - max{X, Y} ' 
In the same way, considering now the ^-coordinate, we get 

Then, considering y fixed, a similar argument gives 

(y + i)(z + i) 



1^1 < 



M - max{X, Y} ' 

□ 



Corollary 14.6. Let B and C he non-zero rational integers with gcd(-B, C) = 1 
and let D be an integer. Define the plane U (with A = 0), S and M as in the above 
Lemma. Then 



M < (X + 1) 1 



Y 
]C\ 



and M < (X + 1) 1 



\B\ 



Moreover, if we suppose that 

M > max{X + Y + 1,X + Z +1} 

then 



\B\ < 



{X + l)iZ+l) 
M-X 



\C\< 



(x + i)(y + i) 

M-X 



Proof. The proof is similar to that of the Lemma, but simpler. We omit the details. 

□ 



Lemma 14.7. Let Ri, Si and Ti be positive integers and consider the set 

ti = {{xi,X2) = {r + sl3i,t + spa) : < r < Ri, < s < Si, <t <Ti} 
and assume that 

Card El = {Ri + l){Si + l)(Ti + 1). 
Let (A, /u) € \ {(0,0)} and let c be a complex number. Let x C' positive real 
number. Then, for any c, the number M of elements {xi,X2) & ^i such that 

Xxi + flX2 = c 

satisfies 

(20) M < max{Ei+5i+l, Si+Ti+l, Ri+Ti+1, x((iii+l)(5i+l)(Ti+l))''''} 

— except if, either there exist two non-zero rational integers ri and si such that 

rib2 = sibi 
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with 

(i?i + l)(Ti + l) 



\ri\< 



and \si I < 



x[iRi + l){Si + l)(ri + l)j - max{i?i,Ti} 

(5i + i)(ri + i) 



x[{Ri + l){Si + l)(ri + l)j - max{5i,ri} 
or there exist rational integers ri, si, ti and t2, with riSi ^ 0, such that 

{tibi + ri 63)51 = ri52t2, gcd(ri,ti) = gcd(si,t2) = 1, 
which also satisfy 

(i?i + l)(5i + 1) 



< |risi| < 6- 
\siti\ <$■ 
and \rit2\ < S ■ 



x[{Ri + l){Si + + l)j - max{i?i, 5i} 

(5i + i)(ri + i) 

x[{Ri + l){Si + l)(Ti + l)j - max{5i, Ti} 

(Ri + i)(ri + 1) 

/ \l/2 

x[{Ri + l){Si + l)(Ti + l)j - max{i?i, Ti} 



where 

S = gcd(ri,si). 

// the previous upper bound ^2(J\) for M holds then, for all (A, /i) G \ {(0, 0)}, we 
have 

CardjAxi +11x2 : {xi,X2) G Si} 

> (i?i + l)(^i + l)(Ti + l) 

~ inax|i?i + 5*1 + 1), 5*1 +T1 + I, i?i +ri + 1, x((i?i + l)(5'i + l)(Ti + 1))^^^} 

Proof. Let 

El = {{r,s,t) e : < r < i?i, < s < 5*1, < i < Ti}. 
Recall the notation 

xi=r + Pis, X2 = t + P3S, /3i = ^, /Jg = ^. 

02 62 

For (A,^) e \ {(0,0)}, we consider the cardinality 

N = CardjAxi + 11x2 : {xi,X2) £ Si}. 

We put 

M = max Card | (xi, X2) G Si : Aa;i + 11x2 — c| 
cec ^ ^ ' 

and 

lie = {(^1, ^2) e : Azi + ^Z2 = c}. 

We clearly have 

TV > Card Si /M, 

so that the last claim of the Lemma is proved and we may also suppose that 1)20(1 
does not hold. 



36 



YANN BUGEAUD, MAURICE MIGNOTTE, SAMIR SIKSEK 



Consider a complex number c such that the number of points {xi,X2) G Si for 
which Xxi + iJbX2 € He is maximal (and so equal to M). 

• If /X = 0: Apply the previous Corollary with (.x, y, z) i— > (r, s, t), {X, Y, Z) i— > 
{Ri,Si,Ti), {A,B,C) ^ (62/^3,61/^3,0), where 

ds = gcd(6i,52), 

and (62/0(3, bi/ds) (ri, si). Then we get the wanted assertion (the 'either' 
case). 

Now we assume /U ^ and, to simplify the notation we take /x = 1. 

• If A = 0: Now, as above, apply the previous Corollary but with {A, B, C) 
(0,63/(^1,62/^1), where 

di = gcd(62,63), 
and (62/(^1,63/^1) I— > (si,t2). Then we get the asserted relation 
[tibi + ri 63)51 = ri62t2 

with ri = 1 and ti = 0, and the asserted bounds on ri, si, t\ and t^,- 

• If A61 + 63 = 0: In this case {A, B, C) ^ (—63/(^2, 0, 61/(^2), where 

d2 = gcd(6i,63), 

and (6i/(i2, —63/6(2) ^ (ri,ti). Then we get the asserted relation 
(ti 61 +ri 63)51 =ri62t2 

with si = 1 and = 0, and the asserted bounds on n, si, ti and t2- 

• If A^(A6i + 63) 7^ 0: Since M > 6*1 + 1, there exist two distinct triples 

(r, SQ,t) and (r', sq, t') G such that 

A(r + /3iSo) + {t + jSsSo) = A(r' + PiSq) + {f + P3S0), 

which gives A(r' — r) = t — t' , where we suppose (as we may) that r is 
minimal (then r' > r) and also that r' — r > is minimal. Put ri = r' — r 
and ti = t — t', then 

A = ti/ri. 

Since M > R-i+1, there exist two distinct triples (ro, s, t) and (ro, s',t') G 
E such that 

^162^0 + (ii6i + ri63)s + ri62i = ti62ro + (ii6i + ^63)5' + ri62t', 
which gives now a relation of the form 

{tihi +ri63)si = ri62t2, 
for which we may suppose that 

gcd(ri,fi) = gcd(si,t2) = 1- 
Now we arc ready to apply the above Corollary with 
{A,B,C) ^ {tis^/5,nt2/5,nsi/5), 

where 

5 = gcd(iisi,rit2,risi), 
and we get the conclusion, except that wc have to prove that 5 = gcd(ri, si). 

Suppose that p is a prime divisor of 6, then p \ nsi. If p | ri then p \ si 
and p] ti, thus p \ riti: contradiction. If p | si then p \ n and p \ ti, thus 
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p \ sit2'- contradiction. Hence, p \ ri and p \ si and p \ t\t2- And now it is 
easy to conclude that 

6 = gcd(ri,si). 
This ends the proof of the Lemma. 

□ 

Remark. Before leaving this Subsection, it is important to notice that the conclu- 
sion of the zero-lemma, namely '. . . the only polynomial P € C[Xi,X2,y] with 
degxi P < for i = 1, 2, and degy P < L which is zero on the set Si + S2 + S3, is 
the zero polynomial' applied to the interpolation matrix considered above implies 
that this interpolation matrix is of maximal rank, which means that there exists a 
determinant A as above which is nonzero. 

14.5. Statement of the main result: a lower bound for the linear form. 

If we gather the results obtained in the previous subsections, we get the following 
theorem. 

Theorem 3. We consider three non-zero algebraic numbers ai, and as, all ^ 1 
which are either all real or all complex of modulus one. Moreover, we assume that 

{either ai, 0:2 and are multiplicatively independent, or 
two multiplicatively independent, the third a root of unity ^ 1 . 

We also consider three non-zero rational integers bi, 62, 63 with gcd(6i, 62, 63) = 1; 
and the linear form 

A = 61 logai + 62 loga2 + bs log 0:3 7^ 0, 

where the logarithms of the ai are arbitrary determinations of the logarithm, but 
which are all real or all purely imaginary. Without loss of generality, we assume 
that 

62I loga2| = \bi logail + I63 log 03] ± |A|. 
Let K, L, R, Ri, R2, R3, S, Si, S2, Ss, T, Ti, T2, T3 be rational integers which 
are all > 3, with 

L>5, R> R1+R2 + R3, S> S1 + S2 + S3, T >Ti + T2+n. 
Let p > 1 be a real number. Assume first that 

(21) + i " ^ ~ H) ^ + 1) + aHaiR + a2S + azT) 

+ V{K-l) log&-21og(e/2), 
where N = K^L, V = [Q{ai,a2,a3) : Q] / [R(ai,a2,a3) : M], e = exp(l), 



(M) 



1 N /^"^ \ ^^^^ 

^=I-T^' b={b2VoKb2Co)[l[ki) 



with 



and 



R-1 , (5- 1)61 ^ T-1 (5-1)63 



ai> p\ logail -log \a,\+2T>h{ai), « = 1> 2, 3. 
If, for some positive real number x, 
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(i) (i?i + l)(^i + l)(Ti + l) > 

+ + Si+T^ + l, R^+T^ + l, x((i?i + l)(^i + l)(Ti + 1)) 

(ii) Card{a>^4 : < r < < s < Si, < i < Ti} > i, 

(iii) (i?2 + l)(5'2 + l)(r2 + l) > 2^2. 

(iv) Card : < r < R2, < s < S2, < t < T2} > 2KL, and 

(v) (i?3 + 1)(53 + 1)(T3 + 1) > 6K^L, 
then either 

A' > p-^^, 

where 

r Li?,e^^l^l/(2fei) L5'e^S|A|/(262) ^ygLT|A|/(2h3) > 

^'='^'•"^^■4 — m — — ' — m — I' 

or at least one of the following conditions (CI), (C2), (C3) hold. 
(CI) |6i| < Ri and I62I < .5i and I63I < Ti, 

(C2) |6i|<i?2 and |62|<5'2 and I63I < T2, 

(C3) either there exist two non-zero rational integers n and si such that 

ri62 = sibi 

with 

(i?i + l)(Ti + l) 



n < 



( \ 1/2 

x((i?i + l)(5i + l)(ri + l)j - niax{i?i,Ti} 

(5i + l)(Ti + l) 



and |si| < 



x((i?i + l)(5i + l)(ri + l)j - max{Si,Ti} 
or there exist rational integers ri, si, ti and t2, with risi 7^ 0, such that 

{tibi + ri 63)31 = rib2t2, gcd(ri, ii) = gcd(si, t2) = 1, 
which also satisfy 

Isiiil < <5. 
and \rit2\ < S ■ 



x((i2i + l)(5i + l)(Ti + l))'^'- 
(5i + l)(Ti + 1) 


niax{i?i, 5i} 


x((i?i + l)(5i + l)(Ti + l))'^'- 
(iii + l)(Ti + 1) 


niax{5i,Ti} 


x((i?i + l)(5i + l)(Ti + l))'^'- 


niax{i?i,Ti} 



where 

6 = gcd(ri,si). 
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Warning . — In the above theorem, the roles of (ai, bi) and (a2, 62) are not com- 
pletely symmetric. Even if we do not make the hypothesis ai > 03 (and, of course, 
do not use it), in practice it is sometimes better to choose the numerotation such 
that fli > as, but one has also to deal with (C3) which is also non-symmetrical. . . 

14.6. An estimate for linear forms in two logarithms. We need to use linear 
forms in two logarithms in a very special situation (related to condition (C3) above) 
and it is difficult to find an easy-to-use result for such a case. This is the reason 
why we write a suitable application of in this Section. We apply the Corollary 
of Theorem 2 of [S]: 

Proposition 14.8. Consider the linear form in two logarithms 

A = &2 log a2 — bi log ai , 

where bi and 62 are positive integers. Suppose that a\ and a2 are multiplicatively 
independent. Put 

I?= [Q(Qi,a2) :Q]/[R(ai,a2) : M]. 
Let ai, 02, h, k be real positive numbers, and p a real number, e'^/^ < P < e"^. Put 
A = log/9, x = h/\ and suppose that x ^ Xo for some number Xo ^ o,nd that 

(22) h > max|7.5,3A,P ^log(^ + ^) + log A + 1.285^ + 0.023| , 



(23) a, > max{4, A,p|logai| -loglofil -H2m(a,)}, (i = l,2), 0102 > 100. 
Put 

i' = 4x + 4+l/x, A = max{ai, 02}. 
Then we have the lower bound 

log|A| > -{Co + ci + C2)iX + h)^aia2, 

where 



i- 



I \ 2x(x + 1) / I 3 y 9 3v \ai 02/ 3v^^/aia2 

and 

A(1.5A + 2h) _ l.llAlog(A(2A + 2hy 



ci = , , s9 > C2 = 



(A + h)^aia2 ' (A -I- h)^aia2 

Proof. The only difference with Theorem 2 of pT| is the definition of the term h. 
Put 

1 f V2 + 2^ / 2(l + xo) , 2A ^ 1 , 1 \ , 4A,/2+^ V 

-^0 •= T 5 + \ n + ^ ~ + ~ + q , «i«2 

A \ 3 y 9 3 \ai 02 J 3y'aia2 I 

and 

(1 + V^^)V^ logx , 3 , ^ 3 , log^ 

fix) = log + — + - + log - + — . 

X — 1 ox[x — 1) 2 4 X — 1 

Then the condition on h in Theorem 2 of jSlj is 

h>v(\og(— + —) + log A + /( [i^ol)) +0-023. 
\ \a2 ai/ J 
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Here we can take xo = 3 and it is easy to check that our present hypotheses imply 
Kq > 195. Since f {x) < 1.285 for x > 195, we get the result. □ 

We notice that ci is a decreasing function of x, for x 1 + we have 

1 

ci < 



2aia2 
We also have 

l.llAlog(aia2(A + ^ i.HA log(aia2(A + /i)^ 



C2 < 



hence 

- e{iX + 'h)'^r ^ e(l + x)(A + I)'/^(ara2)«A0 < ' («^«^)"''" 

(notice that the hypotheses of the above Proposition imply x ^ 3 and X + h >9). 

Using these remarks and simplifying the expression of Co using w > 16 we get 
the simpler estimate. 

Corollary 14.9. With the notation and hypotheses of the above proposition, we 
have the lower bound 



log|A| > -(C^ + c'i + 4)(A + /i)2aia2, 



whe 



1 



4 0.177- (0102) 



-9/10 



1 o ' 

2aia2 

14.7. How to use Theorem |3J Here we assume that condition (M) holds. To 
apply the Theorem, we consider an integer L > 5 and real parameters m > 0, p > 1 
(then one can define the a^) and we put 

K ~ \ mLaia2a3\ , with 771010203 > 2. 

To simplify the presentation, even if we do not really need these conditions, we also 
assume 

TO > 1, and ai, 02, 03 > 1. 

We define 



Ri = 


[010203], 


^1 = 


[C10103J, 


Ti 


= [ciOia2j, 


R2 - 


[020203], 




[C20103J, 


T2 


= [C2aia2j, 


R3 = 


[030203], 


^3 = 


[C30103J, 


T3 


= [C3O1O2J, 



where the parameters ci, C2 and C3 will be chosen so that the conditions (i) up to 
(v) of the Theorem are satisfied. 
Clearly, condition (i) is satisfied if 

(01(010203)^)^^^ > X"^°'i02a3i, Ci • o > 2toL, where o = min{oi, 02, 03}. 
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Condition (ii) is true when 20^010203 • min{ai, 02, 03} > L. Thus, since we suppose 
m and the all > 1, we can take 

ci = max < (xmL)-^''', ( j \. 

To satisfy (iii) and (iv) wc can take 

C2 = max |2^/^(mL)2/^ ^/m/a l| . 

Finally, because of the hypothesis (M), condition (v) holds for 

C3 = [Qm^fl^L. 

Remark. When ai, a2, are multiplicatively independent then it is enough to 
take ci and C3 as above and 

C2 = 2i/3(mL)2/3. 

Then we have to verify the condition (|21(l . When this inequality holds, one obtains 

iA'i > p-^^ 

and we get 

log I A| > -ifL log p - log(max{i?, 5, T) ■ L), 

except maybe if at least one of the conditions (CI), (C2) or (C3) holds. 

Now consider the conditions (CI), (C2) and (C3). For conditions (CI) and 
(C2) we have in particular 

(CI) or (C2) 62 < max{S'i,S'2}. 

Condition (C3) will be studied for an example. Put 

n = 5r[, si = (5s'i, 

where 

5 = gcd(ri,si). 
We just notice, for the second alternative, 

[tibi + ri63)si = rib2t2, gcd(ri,ii) = gcd(si,t2) = 1, 
that r'^ I 5i, say 61 — r[b[, hence 

{tib[+6b3)s[ = b2t2, with bi^rib[. 
If i2 7^ this shows that s[ \ 62, say 62 = s[b2, so that 

tib'i + Sb^ = 62^2, with 61 = r[b'^, and 62 = s'i^2- 

15. Completion of the Proof of Theorem ^ 

Having given our new bounds for linear forms in three logarithms we now use 
them to complete the proof of Theorem ^ We have indeed shown in Lemma 111.11 
that if {x, y,p) is a missing solution then p > po where po is given in Tabled To 
complete the proof it is enough to show that p < po- In Section [T^ we wrote down 
the linear form in logarithms we obtain for each outstanding value of D. We will 
content ourselves by giving the details of this calculation for D = 7. The other 
cases are practically identical (but with different constants, and a different number 
of iterations). 
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We defined 

X - 

A = los 



and we have seen that 

log|A| < ~| logy + log(2.2V7). 

Writing ao = (1 + \/— 7)/2 we saw that the linear form is given by 

A = 21og(ao/ao) +plog(7/7) + 

for some rational integer q (which is not necessarily prime, but we have some lack 
of notation!) with \q\ < p, and we get 

log|A| > -/fLlog/9-log(max{i?,5,T} -L), 

except maybe if at least one of the conditions (CI), (C2) or (C3) holds. 
We have already seen that 

(CI) or (C2) ^ p< max{S'i,S'2}. 

Thus, if p > xnax{Si,S2} then (CI) and (C2) do not hold and then (C3) holds, 
and by the study in the previous section, either 

2(5i + i)(ri + i) 



62 = p < 



X((i?i + l)(5i + l)(ri + 1)) ' - max{5i,ri} 
or we obtain a relation 

t'b' + t"p + g = 0, with 6' = 1 or 2, 
(here we have used the fact that p is prime, which implies |si| = 5 — 1), where 

un< (5i + l)(ri + l) 



x((i?i + l)(5i + l)(ri + 1))'^' - max{5i,Ti} ' 
Hence, 

, ,n kl 2|t'| 1 2\t'\ 
\t"\ < lli + -L-i < - + -^, 
p p 2 p 

which implies 

p< <s. + im + i) 

xiiRi + l){Si + l)(ri + 1))'^' - max{5'i,Ti} 
If the second alternative holds then 

2(^1 + l)(Ti + 1) 



\q\ < 2|t'| < 



X((i?i + l){Si + l)(Ti + 1))'/' - max{5i, Ti} 
and we can apply [25| to the linear form in two logs 

A = Logai +pLoga2, 
taking a2 and Loga2 — log 012 as before, but now 

ai — ±{a/a)'^ and Logai — logcki + iqn. 
And we get for example (using CoroUaire 1 and the notation of 25.) 

log|A| > -31 X logAi X logA2 X (max{21,logp})^. 
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Now we proceed effectively to the computation of an upper bound for p. The 
first step is to recall that we have proved in Lemma fl3.4l by applying Matveev's 
Theorem (Theorem that 

p < 6.81 X 10^^ 

We then apply our Theorem O with the initial condition p < 6.81 x 10^^ and 
with the lower bound 

y>22; 

note that we do not yet assume our lower bound (|14(l obtained through the modular 
approach. There are two reasons for this: 

• The first reason is that we would like to demonstrate how powerful our new 
lower bound for linear forms in three logarithms is, even without the help 
of the modular approach. 

• The second reason is that when we later make the assumption (|14|) . and 
apply our lower bound for linear forms in three logarithms, the reader will 
be able to appreciate the saving brought by the 'modular lower bound' for 

y- 

So for now we assume simply that y > 22 which can be deduced from the fact that 
y is even, is not a power of 2 and that —7 is a quadratic residue for every odd prime 
factor of y (see [IHl)- Applying Theorem 13 we get 

p < 3.05 X 10^ 

with the choices L = 120, p^5,m^ 106.2055121, x = 0.4 and 

Ri = 46385, Si = 54196, Ti = 37763, R2 = 107649, S2 = 125777, T2 = 87639 
and 

i?3 = 765790, S3 = 894748, T3 = 623444 

unless at least one of the conditions (CI), (C2), (C3) holds. For these values, it 
is clear that — since we know that p > 10* — conditions (CI) and (C2) do not 
hold ^. We also see that we must have t" — and then that 

\q\ < 880, 

(which contains also the exceptional case q ~ 0). Using CoroUaire 1 of |25| . we get 
log|A| > -31 logAi \0gA2 X (max{21,logp})^ 

with (here) 

logAi = (|g| +2)tt, \0gA2 = i max{7r,logy}; 

which leads to 

p<8x 10^ 

in contradiction with our hypothesis p > 10*. Thus we have proved that 

p < 3.05 X 10^ 

^To be more precise we can take the above values for Si , ^2 and S3 independently of y but the 
-Ri's and Ti's have to be increased for y > 22, as can be seen on the definition of the parameters 
given in the previous section [ai and as are independent of y but not 02] . Luckily, the larger y 
is, the better our resulting estimate for p will be and thus we can always replace y by some lower 
bound for it. 
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Now we iterate the same process, beginning with this new upper bound on p. After 
four iterates (keeping the choices L = 120 and p = 5), we get 

p < 1.11 X 10^ 

The reader should compare this bound with the bound p < 6.81 x 10""^^ obtained 
by Matveev's Theorem. 

We now assume our 'modular lower' bound for y in (|14|l . and then we obtain the 
much better bound for p (taking L — 115 and p — 5.4) 

p < 3.94 X 10^ 

We try to change a little the linear form. As before we choose e = ±1 such that 
the principal determination of log(e7/7) has an absolute value < 7r/2 and we take 
now s' such that the absolute value of the the principal determination of log(£'a/Q;) 
is minimal and thus < 7r/2. Then we have to distinguish two cases 

(I) &i = 2, ai = e'a/a, 62 = P, ^2 = £7/7, 63 = 9, "3 = -1, 
and 

(II) bi = 2, Qfi = e'a/a, 62 ==9, "2 = -1, ^3 ^ P, "3 = £7/7- 

The study of the first case corresponds exactly to our above study, but with the 
better value 

I log ai I = min | log - , tt - log - \ . 

la a J 

In this case we get now (with the same choices of L, p and x as above) 

p < 1.56 X 10^ 

Concerning case (II), we first notice that 

(CI) or (C2) =^ p< max{Ti,T2}, 

an implication which is essentially equivalent to the previous one whose conclusion 
was p < max{5i, 5*2}. (Indeed, the present T^'s play the role of the previous Si's, 
and both are bounded independently of y.) 

Now we study condition (C3). For the first alternative 

7-162 = S161, 

we get 

< ¥^ [ + 1 ) 

and we can apply |25j to the linear form in two logs 

A= (log(£'a/a)2 + glog(-l)) +plog(e7/7), 

which works quite well. 

Consider now the second alternative, which gives here 

2s't' + r't"q + r's'p = 0. 

The cases t' = and t" = are very easy to treat, we omit the details and assume 
t't" 7^ 0. As before, dividing the above relation by r', we obtain 

s't'b' + t"q + sp 0, with 6' = 1 or 2. 
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We can write 



t'h' + t" q +p = 0, with q = s'q' , 



and 



t'A = g'(log(a7'*"/''') +sYlog(-l)) +plog(47a? 



2/6'^ 



which we can estimate as a hnear form in two logs. Now we have to use Corol- 
lary 114.91 above. 

We have the following data. We choose L ~ 115, p ~ 5.5 and x = 1 and we get 
p < 1.81 X 10^, 



with 



and 



Ri 



when (C3) does not hold, 
117653, Si = 31819, Ti = 19991 



1.03, 



'(5i + i)(ri + i) 



(i?i 



1) 



= 76 



^2 



1.03, 



'(i?i + i)(ri + i) 



= 276. 



(^1 + 1) 

we find p < 9 x 10'' when (C3) holds. Finally, 



Using Corollary 14.8 for p = 
we have proved that 

p < 1.81 X 10^ 

Notice that here the use of the Corollary 1 of [21] produces a result which is too 
weak for our purpose, this is the reason why we have written here a special lower 
bound for linear forms in two logarithms. 

Remark. We notice that without the modular lower bound for y we were able to 
show that p < 1.11 X 10^ whilst with this modular lower bound we were able to 
improve this to p < 1.81 x 10^. Whilst it is certainly possible to reach the former 
target with the methods of this paper, it would have taken about 6 times as long 
as it took to reach the latter. From this it is a plausible guess that without the 
modular lower bound for y the computational part for the entire proof for all the 
values ol 2 < D < 100 might have taken at least 1200 days rather than 206 days. 
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16. Tables 



L) 


Solutions , |?/| , Ti) 


1 


(0.1, ;0 


2 


(5,3,3) 


3 




4 


(2, 2, 3), (11, 5, 3) 


5 




6 




7 


(1, 2, 3), (181, 32, 3), (3, 2, 4), (5, 2, 5), (181, 8, 5), (11, 2, 7), (181, 2, 15) 


8 


(0,2,3) 


9 




10 




11 


(4, 3, 3), (58, 15, 3) 


12 


(2,2,4) 


13 


(70,17,3) 


14 




15 


(7, 4, 3), (1,2, 4), (7, 2, 6) 


16 


(0.2,4). (4,2.5) 


17 


(8.3,4) 


18 


(3,3,3), (15,3,5) 


19 


(18,7,3), (22434, 55,5) 


20 


(14,6,3) 


21 




22 




23 


(2, 3, 3), (3, 2, 5), (45, 2, 11) 


24 




25 


(10,5,3) 


26 


(1,3,3), (207, 35,3) 


27 


(0,3,3) 


28 


(6, 4, 3), (22, 8, 3), (225, 37, 3), (2, 2, 5), (6, 2, 6), (10, 2, 7), (22, 2, 9), (362, 2, 17) 


29 




30 




31 


(15,4,4), (1,2,5), (15,2,8) 


32 


(7, 3, 4), (0,2, 5), (88, 6, 5) 


33 




34 




35 


(36,11,3) 


36 




37 




38 




39 


(5, 4, 3), (31, 10, 3), (103, 22, 3), (5, 2, 6) 


/in 


('CO 1 4 Q"\ 


41 




42 




43 




44 


(9,5,3) 


45 


(96, 21, 3), (6, 3, 4) 


46 




47 


(13, 6, 3), (41, 12, 3), (500, 63, 3), (14, 3, 5), (9, 2, 7) 


48 


(4, 4, 3), (148, 28, 3), (4, 2, 6) 


49 


(524, 65, 3), (24, 5, 4) 


50 
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D 


Solutions (Ixl, Ivl, n) 


51 




52 




53 


(26, 9, 3), (156, 29, 3), (26, 3, 6) 


54 


fl7 7 3) 


55 


f3 4 3) f419 56 3),f3,2,6) 


56 


f76, 18,3), (5,3,4) 


57 




58 




59 




60 


f2 4,3) (1586,136 3) (14 4 4), (50354 76 5), (2, 2, 6), (14, 2, 8) 


61 


(8,5,3) 


62 




63 


(1.4, 3). (13537, 568. 3), (31. 4. 5), (1. 2. 6), (31, 2, 10) 


64 


(0.4,3). (0,2.0), (8.2,7) 


65 


(4,3,4) 


66 




67 


(110,23,3) 


68 




69 




70 




71 


(21 8 3) (35 6 4), (46,3,7), (21,2,9) 


72 


(12,6,3), (3,3,4) 


73 




74 


(985,99,3), (13,3,5) 


75 




76 


(7,5,3), (1015, 101,3) 


77 


(2,3,4) 


78 




79 


(89,20,3), (7,2,7) 


80 


(1,3,4) 


81 


(46, 13,3), (0,3,4) 


82 




83 


(140,27, 3) 


84 




85 




86 




87 


(16, 7, 3), (13, 4, 4), (13, 2, 8) 


88 




89 


(6 5,3) 


90 




91 




92 


(6, 2, 7), (90, 2, 13) 


93 




94 




95 


(11, 6, 3), (529, 6, 7) 


96 


(23,5,4) 


97 


(48,7,4) 


98 




99 


(12.3,5) 


100 


(5, 5, 3), (30, 10, 3), (198, 34, 3), (55, 5, 5) 
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